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Abstract. We exhibit a probabilistic symbolic algorithm for solving zero- 
dimensional sparse systems. Our algorithm combines a symbolic homotopy 
procedure, based on a flat deformation of a certain morphism of affine varieties, 
with the polyhedral deformation of Huber and Sturmfels. The complexity of 
our algorithm is quadratic in the size of the combinatorial structure of the 
input system. This size is mainly represented by the mixed volume of Newton 
polytopes of the input polynomials and an arithmetic analogue of the mixed 
volume associated to the deformations under consideration. 



1. Introduction 

Numeric and symbolic methods for computing all solutions of a given zero- 
dimensional polynomial system usually rely on deformation techniques, based on a 
perturbation of the original system and a subsequent (numeric or symbolic) path- 
following method (see, e.g., PP, 0, PS])- The complexity of such 
algorithms is usually determined by geometric invariants associated to the family 
of systems under consideration (see, e.g., [E], [2FJ, [S|, El, E2, EH, [El, [E3j, 
31 , 39 ), typically in the form of a suitable (arithmetic or geometric) Bezout 
number (see 051, [2E], [HSj, ESI, [23, USD- 
Sparse elimination theory is concerned with finding bounds for such Bezout 
numbers in the case of a sparse polynomial system. Its origins can be traced 
back to the results by D.N. Bernstein, A.G. Kushnirenko and A.G. Khovanski (E], 
|29j . |28j ) that bound the number of solutions of a polynomial system in terms 
of certain combinatorial invariants. More precisely, the Bernstein-Kushnirenko- 
Khovanski (BKK for short) theorem asserts that the number of isolated solutions 
in the n-dimensional complex torus (C*) n of a polynomial system of n equations 
in ri unknowns is bounded by the mixed volume of the family of Newton polytopes 
of the corresponding polynomials. 

Numeric (homotopy continuation) methods for sparse systems are typically based 
on a family of deformations called polyhedral homotopies ([25J, [SI], [S3]). Poly- 
hedral homotopies preserve the Newton polytope of the input polynomials and 
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yield an effective version of the BKK theorem (see e.g. [23, |26|~l. More pre- 
cisely, suppose that we are given a zero-dimensional (Ai, . . . , A n )-sparse system 
defined by n polynomials ft, . . . , f n in n variables, where Ai, . . . , A„ are the sup- 
ports of ft, ...,/ n , and let V C (C*) n be the variety defined by the common 
zeros of ft,...,f n over (C*)™. Then a polyhedral homotopy consists in an alge- 
braic curve W C (C*) n+1 such that the projection tt : W — * C* onto the first 
coordinate is dominant with generically finite fibers whose degree is the mixed vol- 
ume MV(conv(At), . . . , cont>(A n )) of the convex hulls of Ai, . . . , A„, the identity 
7r (1) = {1} x V holds and the first terms of the Puiseux expansions of the branches 
of W lying above can be easily computed. Numerical continuation methods com- 
pute the first terms of these Puiseux expansions and then follow the branches of 
W along the interval [0, 1] to obtain approximations to all the points of the input 
variety V. 

From the symbolic point of view, a family of homotopy algorithms is based on 
a flat deformation of a certain morphism of affine varieties. This deformation, 
implicitly considered in the papers (20j, is isolated in [21] and refined in |5()j . 

0, in order to solve particular instances of a parametric system with 
a finite generically unramified linear projection of low degree. More precisely, let 

V be a zero-dimensional variety of C™, defined by a "square" system f\ = ■ ■ ■ = 
f n = 0, and let be given an algebraic curve W C C ,l+1 and a (dominant) projection 
mapping tt : W — > C which represent a deformation of V. Then, from a complete 
description of a generic fiber of the projection tt : W — ► C, it is possible to compute 
a complete description of the input fiber, say 7r -1 (l) = {1} x V. The complexity of 
this procedure can be roughly estimated by the product of two geometric invariants: 
the degree of the morphism tt and the degree of the curve W. The algorithm is 
nearly optimal in worst case and has good performance over certain well-posed 
families of polynomial systems of practical interest (see [21], [10], 0, |12|'l. 

In this article we combine these symbolic techniques, particularly in the version of 
[7j, with the polyhedral deformation |25| . in order to derive a symbolic probabilistic 
algorithm for solving sparse zero-dimensional polynomial systems with quadratic 
complexity in the size of the combinatorial structure of the input system. More 
precisely, suppose that we are given polynomials ft, ■ ■ ■ , f n of Q[Xt, ■ ■ ■ , X n ] such 
that the system ft — 0, . . . , f n — defines a zero-dimensional affine subvariety 

V of C n . Denote by At,..., A„ C Z" the supports of ft, ■ ■ ■ , f n , an d assume 
that G Ai for 1 < i < n and the mixed volume D of the Newton polytopes 
conv(At), . . ■ , conv(A n ) is nonzero. Then, given a "sufficiently generic" lifting 
function u) for (Ai, . . . , A„) and the corresponding mixed subdivision, we exhibit 
an algorithm which computes a complete description of the solution set V of the 
input system ft = 0, . . . , f n = 0. 

The polyhedral deformation under consideration requires that the coefficients of 
the input polynomials satisfy certain generic conditions (see Section l3.1|) . For this 
reason, we introduce some auxiliary generic polynomials g\ , . . . , g n with the same 
supports Ai, . . . , A„ and consider the perturbed polynomial system hi := ft + <?; 
for 1 < i < n. We first solve the system ht = 0, . . . , h n = and then recover the 
solutions to the input system ft = 0, ...,/„ = by considering a homotopy of type 
ft + (1 - T)gt, ...,/„ + (1 - T)g n (in Section©. 

The system ht = 0, . . . , h n = is solved by considering the polyhedral homotopy 
of |25| . This homotopy introduces a new variable T and deforms the polynomial 
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hi by multiplying each nonzero monomial of hi by the power of T determined by 
the given lifting function uj. From the gcncricity of the coefficients of hi, . . . , h n we 
conclude that the roots of the resulting parametric system hi = 0, . . . , h n = are 
algebraic functions of the parameter T whose expansions as Puiseux series can be 
obtained by "lifting" the solutions to certain zero-dimensional polynomial systems 
fvQ = ■ ■ • = /irv7 = associated to the lower facets C 7 of the lifted polytopes 
conu(Ai), . . . , conv(A n ) defined by u). These polynomial systems hfQ = ■■■ = 

h^l = can be easily solved due to their specific structure, which enables us to use 
their solutions as a starting point for our computations (see Section l4~Tl for details). 

The complexity of our algorithm is mainly expressed in terms of two quanti- 
ties related to the combinatorial structure of the input system: the mixed volume 
D := M (conv(Ai), . . . , conv(A n )} and certain (nonarchimedean) heights E, E' as- 
sociated to our polyhedral deformations. These heights, which are an arithmetic 
analogue of the mixed volume D (see 0T], 02), can be bounded in terms of certain 
mixed volumes associated to the polyhedral deformation under consideration, with 
equality for a generic choice of the coefficients of the polynomials hi (see Lemma l2"3l 
below; compare also with [431 Theorem 1.1]). Therefore, we may paraphrase our 
complexity estimate as saying that it is quadratic in the combinatorial structure of 
the input system, with a geometric and an arithmetic component. More precisely, 
our algorithm requires Cn°^D max{i?, E'} arithmetic operations over Q (up to 
polylogarithmic terms), where C is the number of arithmetic operations required 
to evaluate the polynomials h% and /, + g%, and E and E' denote the height of the 
varieties defined by h x =0,...,h n = and A + (1-T) 5l = 0, . . . , f n +(l-T)g n = 
respectively. 

This improves and refines the estimate of |7] in the case of a sparse system, 
which is expressed as a fourth power of D and the maximum of the degrees of the 
varieties hi = 0, . . . ,h n = and fi + (1 - T)gi = 0, ...,/„ + (1 - T)g n = 0. We 
observe that this maximum is an upper bound for the heights E and E' respectively. 
On the other hand, it also improves 01], 0S]> which solve a sparse system with a 
complexity which is roughly quartic in the size of the combinatorial structure of the 
input system. Finally, we provide an explicit estimate of the error probability of 
all the steps of our algorithm. This might be seen as a further contribution to the 
symbolic stage of the probabilistic seminumeric method of [23], which lacks such 
analysis of the error probability. 



2. Notions and notations 

2.1. Sparse Elimination. Here we introduce some notions and notations of con- 
vex geometry and sparse elimination theory (see e.g. |18j . [25| ) that will be used in 
the sequel. 

Let Xi, . . . ,X n be indeterminates over Q and write X := (Xi, . . . ,X n ). For 
q:=(qi,..., q n ) G Z n , we use the notation X q := Xf ■ ■ ■ X%> . Let / := J2 q c q x<! 
be a Laurent polynomial in Q[X, X^ 1 ] := Q[X\, Xf , . . . , X n , X~ x ]. By the support 
of / we understand the subset of Z" defined by the elements q G Z™ for which c q ^ 
holds. The Newton polytope of / is the convex hull of the support of / in E™. 
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A sparse polynomial system with respect to a priori fixed finite subsets Ai, . . . , A„ 
of (Z>o)" is defined by polynomials 

fi(X) := ^s xq (l<*<n), 

qeAi 

with ai iQ £ C \ {0} for each q £ Ai and 1 < i < n. 

For a finite subset A of Z™, we denote by Q := conti(A) its convex hull in K n . 
The usual Euclidean volume of a polytope Q in K™ will be denoted by vol^(Q). 

Let Qi, . . . , Q n be convex polytopes in R™. For Ai, . . . , A n <E R>o, we use the 
notation XxQx + • • • + X„Q n to refer to the Minkowski sum AiQi + • • • + X n Q n '■= 
{x 6 K" : 3xi € Qi, . . • , 3x n € Q n such that x = \\Xi + • • • + X„x n }. Consider 
the real-valued function (Ai, . . . , A„) i— > woZRn(AiQi + • • • + X n Q n ). This is a ho- 
mogeneous polynomial function of degree n in the Ai (see e.g. Chapter 7, 
Proposition §4.4.9]). The mixed volume MV(Qx, . . . , Q n ) of Qx, . . . ,Q n is defined 
as the coefficient of the monomial Ai - •• A n in volRn(XxQx + • • • + X n Q n ). 

For i = 1, . . . ,n, let A, be a finite subset of Z™ and let Qi := conv(Ai) de- 
note the corresponding polytope. Let /i,...,/ n be a sparse polynomial system 
with respect to Ai,...,A n . The BKK theorem (gj, [213, EE]) asserts that the 
system fx = 0, ...,/„ = has at most MV(Qi, . . . , Q n ) isolated common solutions 
in the n-dimensional torus (C*) ra , with equality for generic choices of the coeffi- 
cients of fx, ...,f n . Furthermore, if the condition e Qi holds for 1 < i < n, 
then MV(Qx, ■ ■ ■ , Q n ) bounds the number of solutions in the n-dimensional affinc 
complex space A" := A n (C) (see [33]). 

Assume that the union of the sets Ai, . . . , A„ affinely generate Z™, and consider 
the partition of Ai, . . . , A„ defined by the relation Aj ~ Aj if and only if A,; = Aj. 
Let sgN denote the number of classes in this partition, and let A^ 1 ', . . . , A^ C Z™ 
denote a member in each class. Write A := (Jv- 1 ', • • ■ , A^). For t = 1, . . . , s, 
let kg :— #{j : A^ = A^}. Without loss of generality, we will assume that 
Ax = ■ ■ ■ = A kl = A {1 \ A fel+ i = ■ • • = A fcl+fc2 = A& and so on. 

A cell of A is a tuple C = (C*W, . . . , C^) with CW ^ and C A {1) for 
f < I < s. We define 

type(C) := (dim(com;(C (1) )), . . . , dim(conv(C {s) ))), 

conv(C) :=conv(C w +■■■ + C (s) ), 

#(C) :=#(£«) + + 

i>o?Rn(C) :— vol^ (conv{C)). 

A face of a cell C is a cell C = (C«, . . . ,C^ s) ) of C with C for 1 < I < s 
such that there exists a linear functional 7 : M" — > M that takes its minimum over 
at for f < € < s. One such functional 7 is called an inner normal of C. 
A mixed subdivision of A is a collection of cells £ = {Cx, . . . , C m } of „4 satisfying 
conditions (l)-(4) below: 

(f) cYvax(conv{C j)) = n for 1 < j < m, 

(2) the intersection conv(d) (~l conv(Cj) C K n is either the empty set or a face 
of both conv(Ci) and conv(Cj) for 1 < £ < J < m, 

(3) U^Li conv{Cj) = conv(A), 

(4) 21=1 dim(conw(C^)) — n for f < j < m. 
If £ also satisfies the condition 
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(5) #(Cj) = n + s for 1 < j < m, 
we say that C is a fine-mixed subdivision of A. Observe that, as a consequence of 



conditions (4) and (5), for each cell C 3 = (cf \ . . . , C) s> ) in a fine-mixed subdivi- 



sion the identity dim(conv(cf } )) = #Cf> - 1 holds for 1 < £ < s. 

We point out that a mixed subdivision £ of A enables us to compute the mixed 
volume of the family Q\ = conv(Ai), . . . , Q n = conv(A n ) by means of the following 
identity (see |2S1 Theorem 2.4.]): 



-,(0 



(2.1) 



MV(Qi 



i Qn) 



E 

type(C i ) = (fci,...,fc 3 ) 



fell. .. k s \ ■ VOlRn(Ci). 



A fine-mixed subdivision of A can be obtained by means of a lifting process 
as explained in what follows. For 1 < I < s, let tug : A™> — » R be an arbitrary 
function. The tuple u> := (u>\, . . . ,uj s ) is called a lifting function for A. Once a 
lifting function uj is fixed, the graph of any subset of A^ will be denoted 
by := {(q,u t (q)) £ W l+1 : q E C w }. Then, for a sufficiently generic lifting 
function cj, the set of cells C of A satisfying the conditions: 

(i) dim(com;(CW + • • • + CM)) = n, 

(ii) (C^ , . . . , C^) is a face of (A*- 1 -*, . . . , A^) whose inner normal has positive 
last coordinate, 

is a fine-mixed subdivision of A (see Section 2]). More precisely, we have the 
following result (cf. [23 Section 2]): 

Lemma 2.1. The lifting process associated to a lifting function u> yields a fine- 
mixed subdivision of A if the following condition holds: for every r 1; . . . ,r s £ Z> 
with X)!=i rg > n and every cell (C^\ . . . , C*- 5 -*) with C^' := {qe,o, ■ ■ ■ , qt,r t 
(l<i< s), if 

( ?i,i-gi,o \ / 91,1—91,0 ^1(91,1) — wi(gi,o) \ 



V(C) 



qx,n — 9i,o 



9«, l — 9s, o 



V 9s 



and V{C) 



then 



9s,o / 



rank(V(C)) 



9l,ri ~ 91,0 ^l(9l,ri) - ^l(9l,o) 



q s>1 - q a<0 v s (q s ,i) - ^s(9s,o) 



\ q s ,r s -q s ,o Us(q a ,rs) - u s (q St0 ) I 



rank(y(C)) =n+l. 



Proof. Notice that (l)-(3) are automatically satisfied by the set of cells defined by 
conditions (i)—(ii). Assume that the condition of the statement of the lemma is 
met and consider a cell C — (C^ , . . . , C^) of A satisfying conditions (?) and (ii) 
above. Being C a lower facet of A, the identity dim(con?;(C' 1 ' + • • ■ + C*W)) = 
dim(conv{C^ + ■ ■ • + <5( s >)) must hold. Write = {q efil qe, re } for 1 < £ < s. 
Then we have that rank(V(C)) — dim(< qij — qg^ : 1 < I < ri, 1 < j < r$ >) = 
dim(conv(C^+--- + C^)) = n and rank(V(C)) = dim(conv(C (1 ^ +■ ■ - + C^)) = 
n. Now, the condition stated on uj implies that Yle=i r f- — n an d, taking into 
account that the J2e=i r £ manv vectors qej — qe.o (1 < £ < s, 1 < j < re) span 
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a linear space of dimension n, we conclude that the equality r ^ = n holds, 

which shows that condition (5) in the definition of a fine-mixed subdivision is met. 
Moreover, as dim (com; (C 1 ^)) > dim{conv(C^' + ■ ■ ■ + C*- 5 ' 1 )) for an arbitrary 

cell C, we see that dim(conv(C^)) — rg holds for every 1 < I < s, which implies 
that condition (4) is also valid. This finishes the proof of the lemma. □ 

Note that the condition rank(V(C)) = n+1 can be restated as the non-vanishing 
of the maximal minors of the matrix V(C). Since rank(l^(C)) = n, these maximal 
minors are nonzero linear forms in the unknown values LOg(qi t j) of the lifting func- 
tion. Thus, if Afi = #A {1) for every 1 < i < s, a sufficiently generic lifting function 
can be obtained by randomly choosing the values ujg(qg^) of w at the points of Jv® 
from the set {1,2,..., p2 J ^ 1 H ^ } , with probability of success at least 1 — 1/ p for 

peN. 

In the sequel, we shall assume that a sufficiently generic lifting function and the 
induced fine-mixed subdivision of A are given. 

2.2. Complexity model and complexity estimates. In this section we describe 
our computational model and briefly mention efficient algorithms for some basic 
specific algebraic tasks. 

2.2.1. Complexity model. Algorithms in computational algebraic geometry are usu- 
ally described using the standard dense (or sparse) complexity model, i.e. encod- 
ing multivariate polynomials by means of the vector of all (or of all nonzero) co- 
efficients. Taking into account that a generic n-variate polynomial of degree d 
has ( d ^ n ) = 0(d n ) nonzero coefficients, we see that the dense representation of 
multivariate polynomials requires an exponential size, and their manipulation usu- 
ally requires an exponential number of arithmetic operations with respect to the 
parameters d and n. In order to avoid this exponential behavior, we are going 
to use an alternative encoding of input, output and intermediate results of our 
computations by means of straight -line programs (cf. (23 7 |52j . |38|. [lip. A 
straight-line program j3 in Q(X) := Q(Xi,...,X n ) is a finite sequence of ratio- 
nal functions (/1, ■ • • , /fe) E Q(X) k such that for 1 < i < k, the function fi is an 
element of the set {X\, . . . ,X n }, or an element of Q (a parameter), or there ex- 
ist 1 <i\, i% <i such that fi = fi t f i2 holds, where Oj is one of the arithmetic 
operations +, — , X,-j-. The straight-line program is called division-free if Oj is 
different from 4- for 1 < i < k. A natural measure of the complexity of (3 is its time 
or length (cf. 8 , 48 ), which is the total number of arithmetic operations per- 
formed during the evaluation process defined by (3. We say that the straight-line 
program j3 computes or represents a subset S of Q(X) if S C ■■■,/&} holds. 

Our model of computation is based on the concept of straight-line programs. 
However, a model of computation consisting only of straight-line programs is not 
expressive enough for our purposes. Therefore we allow our model to include de- 
cisions and selections (subject to previous decisions). For this reason we shall also 
consider computation trees, which are straight-line programs with branchings. Time 
of the evaluation of a given computation tree is defined similarly to the case of 
straight-line programs (see e.g. [S3], for more details on the notion of compu- 
tation trees). 
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2.2.2. Probabilistic identity testing. A difficult point in the manipulation of multi- 
variate polynomials given by straight -line programs is the so-called identity testing 
problem: given two elements / and g of C[X] :— C[Xi, . . . , X n ], decide whether / 
and g represent the same polynomial function on C". Indeed, all known determinis- 
tic algorithms solving this problem have complexity at least maxjdeg/, degg}^ 1 ). 
In this article we are going to use probabilistic algorithms to solve the identity 
testing problem, based on the following result: 

Theorem 2.2 ([2], @H|)- Let f be a nonzero polynomial of C[X] of degree at 
most d and let S be a finite subset ofC Then the number of zeros of f in S n is at 
most 

For the analysis of our algorithms, we shall interpret the statement of Theo- 
rem |^21 in terms of probabilities. More precisely, given a fix nonzero polynomial / 
in C[Xl, . . . , X n ] of degree at most d, we conclude from Theorem l2.2l that the prob- 
ability of choosing randomly a point a £ S n such that f(a) =0 holds is bounded 
from above by d/#S (assuming a uniform distribution of probability on the ele- 
ments of 5"). 

2.2.3. Basic complexity estimates. In order to estimate the complexity of our pro- 
cedures we shall frequently use the notation M(m) := mlog 2 m log log m. Here and 
in the sequel log will denote logarithm in base 2. Let R be a commutative ring of 
characteristic zero with unity. We recall that the number of arithmetic operations 
in R necessary to compute the multiplication or division with remainder of two uni- 
variate polynomials in R[T] of degree at most m is 0(M(m)/log(m)J (cf. 
Multipoint evaluation and interpolation of univariate polynomials of R[T] of degree 
m at invertible points ai, . . . , a m £ R can be performed with 0(M(m)) arithmetic 
operations in R (see e.g. [§]). 

If R = k is a field, then we shall use algorithms based on the Extended Eu- 
clidean Algorithm (EEA for short) in order to compute the gcd or resultant of 
two univariate polynomials in k[T] of degree at most m with 0(M(m)) arithmetic 
operations in k (cf. [S]). We use Pade approximation in order to compute 

the dense representation of the numerator and denominator of a rational function 
/ = p/q € with maxjdegp, degg} < m from its Taylor series expansion up to 
order 2m. This also requires 0(M(m)) arithmetic operations in k ([23, [HI)- 

For brevity, we will denote by f2 the exponent that appears in the complexity 
estimate 0(n ) for the multiplication of two (n x n)-matrices with coefficients in 
Q. We remark that the (theoretical) bound < 2.376 is typically impractical and 
we prefer to take f2 := log 7 ~ 2.81 (cf. 0). 

2.3. Geometric solutions. The notion of a geometric solution of an algebraic 
variety was first introduced in the works of Kronecker and Konig in the last years 
of the XlXth century. Nowadays, geometric solutions are widely used in computer 
algebra as a suitable representation of algebraic varieties, especially in the zero- 
dimensional case. 

Let K denote an algebraic closure of a field K of characteristic zero, let A n (K) 
be the n-dimensional space 'R n endowed with its Zariski topology, and let V — 
■ ■ ■ , C } be a zero-dimensional subvariety of A n (K) defined over K. A 
geometric solution of V consists of 
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• a linear form u = u\ X\ + • • • + u n X n G K[X] which separates the points of 
V, i.e. satisfying tt(£W) ^ u(^) if i ^ k, 

• the minimal polynomial m u := IIk^dC^ ~ u (£,^)) £ ^\Y] of u in V 
(where Y is a new variable), 

• polynomials W\ , . . . , w n G -K"|X] with deg Wj < D for every 1 < j < n 
satisfying 

V = {(w 1 (ri),...,w n (Ti)) e iT /rjeK, m u {rj) = 0}. 

In the sequel, we shall be given a polynomial system fx = ■ ■ ■ = /„ = of n- 
variate polynomials of Q[X] defining a zero-dimensional affine variety V G A" := 
A n (C). We shall consider the system fx = ••• = /„ = (symbolically) "solved" if 
we obtain a geometric solution of V as defined above. 

This notion of geometric solution can be extended to equidimensional varieties 
of positive dimension. For our purposes, it will be sufficient to consider the case of 
an algebraic curve defined over Q. 

Suppose that we are given a curve V C A™ +1 defined by polynomials fx, ■ ■ ■ , f n £ 
Q[X, T]. Assume that for each irreducible component C of V, the identity 1(C) fl 
Q[T] = {0} holds. Let u be a nonzero linear form of Q[X] and ir u : V — > A 2 
the morphism defined by Tr u (x,i) := (t,u(x)). Our assumptions on V imply that 
the Zariski closure n u (V) of the image of V under ir u is a hypersurface of A 2 
defined over Q. Let Y be a new indeterminate. Then there exists a unique (up 
to scaling by nonzero elements of Q) polynomial M u £ Q[T, Y] of minimal degree 
defining ir u (V). Let m u £ Q(T)[Y] denote the (unique) monic multiple of M u with 
degy(m ll ) = degy(M„). We call m u the minimal polynomial of u in V. In these 
terms, a geometric solution of the curve V consists of 

• a generic linear form u G Q[X], 

• the minimal polynomial m u G Q(T)[y], 

• elements v u ...,v n of Q(T)[Y] such that ^X { = v t mod Q(T) ® Q[V] 
and degy(wi) < deg Y (m u ) holds for 1 < i < n. 

We observe that degy mu equals the cardinality of the zero-dimensional variety 
defined by fx, . . . , /„ over A" (Q(T)) ■ 

In the sequel, we shall deal with curves V := V(fx, ■ ■ ■ , f n ) C A n+1 as above. 
The complexity of the algorithms for solving the systems /i = ••• = /«, = 
defining such curves will be expressed mainly by means of two discrete invariants: 
the degree and the height of the projection ir : V — > A 1 . The degree of ir is defined 
as the degree deg m u = degy M u of the minimal polynomial of a generic linear form 
u G Q[-Xi, . . . , X n ] and can be considered as a measure of the "complexity" of the 
curve V. On the other hand, the height of ir is defined as deg T M u and may be 
considered as a measure of the "complexity of the description" of the curve V . 

In the sparse setting, we can estimate degy M u and deg T M u in terms of com- 
binatorial quantities (namely, mixed volumes) associated to the polynomial system 
under consideration (see also |4Hj ) . 

Lemma 2.3. Let assumptions and notations be as above. For 1 < i < n, let 
Qi C K" be the Newton polytope of fi, considering fi as an element of Q(T)[X]. 
Let Qx, ■ ■ ■ , Q n G M" +1 be the Newton polytopes of fx, ■ ■ ■ , fn; considering fx, . . . , f n 
as elements of Q[X,T], and let A G be the standard unitary simplex in the 

plane {T = 0} ; i.e., the Newton polytope of a generic linear form u G Q[A"]. Assume 
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that G Qi for every 1 < i < n. Then the following estimates hold: 

(2.2) deg Y M u <MV n (Qx,...,Q n ), deg T M u < MV n (A, Q u . . . , Q n ). 

Furthermore, if there exist c\, . . . , c n G M>o such that Qi C Qi x [0, Cj] for 1 < i < n, 
i/ien t/ie following inequality holds: 

n 

(2.3) deg T M u <^c i AfV , (A,Qi,... ) Q i _i ) Qi +1 ,...,Q n ). 

i=l 

Proof. The upper bound deg F M u < M V^(Qi, . . . , Q n ) follows straightforwardly 
from the BKK bound and the afhne root count in |33|. 

In order to obtain an upper bound for deg T M„, we observe that substituting 
a generic value y G Q for Y we have deg T M U (T,Y) = deg T M u (T,y) = #{t G 
C; M M (t, y) = 0}. Moreover, it follows that M M (t, y) = if and only if there exists a 
point x G A™ with y = it(x) and (x, t) G V. Thus, it suffices to estimate the number 
of points (x, t) G A n+1 satisfying u{x) — y = 0, fi(x,t) — 0, . . . , f n {x, t) = 0. Being 
u a generic linear form, the system 

(2.4) u(X) - y = 0, fx(X, T) = 0, . . . , f n (X, T) = 

has finitely many common zeros in A™ + . Combining the BKK bound with the 
affine root count of jSH] w e see that there are at most MV(A, Q%, . . . , Q n ) solutions 
of l|2.4|l . We conclude that deg T Af u < MV(A, Qi, . . . , Q n ) holds, which shows 

In order to prove l|2.3|l . we make use of basic properties of the mixed volume 
(see, for instance, |171 Ch. IV]). Since Qi C Qi x [0, cj holds for 1 < i < n, by the 
monotonicity of the mixed volume we have 

MV(A, Q u . . . , Q n ) < MV(A, Q 1 x [0, ci], ...,Q n x [0, c n ]). 

Note that Qi x [0, a] = + Si t i, where S^o = Qi x {0} and S^i = {0} x [0, Cj] 
for i = 1, . . . , n. Hence, by multilinearity, 

(2.5) MV(A,Q 1 x[0,ci},...,Q n x[0,c n ])= J2 MV{A, Si tjl , . . . , S n>jn ). 

Cii,...,j'»)e{o,i} n 

If the vector (ji, . . . , j n ) has at least two nonzero coordinates, then two of the sets 
S\_j 1 , . . . , S n j n are parallel line segments; therefore, MV(A, Si i j 1 , . . . , S n ,j n ) = 0. 
On the other hand, if ji is the only nonzero coordinate, the corresponding term in 
the sum of the right-hand side of l|2.5|) is 

MV n+1 (A, Q x x {0}, . . . , Q^ x {0}, {0} x [0, a], Q w x {0}, . . . , Q n x {0}) 

= Ci MV n (A, Qx,.-., Qi-i, Qi+i, • • ■ , Qn)- 

Finally, for (jx j n ) = (0, . . . , 0) we have MV n+1 {A, Q x x {0}, ...,Q n x {0}) = 
since all the polytopes are included in an n-dimcnsional subspace. 

We conclude that the right-hand side of Ij2.5|l equals the right-hand side of H2.3JI . 
This finishes the proof of the lemma. □ 

From the algorithmic point of view, the crucial step towards the computation of 
a geometric solution of the variety V consists in the computation of the minimal 
polynomial m u of a generic linear form u which separates the points of V . In the 
remaining part of this section we shall show how we can derive an algorithm for 
computing the entire geometric solution of a zero-dimensional variety V defined 
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over Q from a procedure for computing the minimal polynomial of a generic linear 
form u (cf. 0, EH, |SD|)- 

Let A := (Ai, . . . , A n ) be a vector of new indeterminates and let K := Q(A). 
Denote by Ik the ideal in K\X\, . . . , X n ] which is the extension of the ideal / := 
I(V) C Q[Xi, . . . , X n ] of the zero-dimensional variety V, and denote by B := 
K[X\, . . . , X n ]/lK the corresponding zero-dimensional quotient algebra. Write 

v !i : 

Set U := AiXi + ■■■ + A n X n e K[X X , ...,X n ] and let m c/ (A, Y) = ]jf=i ( Y - 
U(£(^)) G Q[A, Y] be the minimal polynomial of the linear form U in the extension 
K B. Note that degmjj = D holds, and that dmjj/dY is not a zero divisor 
in Q[A" x V]. Furthermore, mu(A,U) e I(A n x V) C Q[A, X u . . . ,X n ] holds. 
Since I(A n x V) is generated by polynomials in Q[X\, . . . ,X n ], taking the partial 
derivative of mj/(A, U) with respect to the variable A& for 1 < k < n, we conclude 
that 

(2-6) ^(A,C/)X fe + ^(A,C/ )e J(A«xn 

Observe that the degree estimate deg y (dmu / dAk) < D — 1 holds. 

Assume that a linear form u — u\X\ + • ■ • + u n X n £ Q[Xi, . . . , X n ] which 
separates the points of V is given. Substituting uu for A^ in the polynomial 
mjj{A,Y) we obtain the minimal polynomial m u (Y) of u. Furthermore, making 
the same substitution in the polynomials {dmjj / dY){A, Y) Xk + (dmu/dAk)(A, Y) 
of (|2.(ill for 1 < k < n and reducing modulo m u (Y), we obtain polynomials 
(dm u /dY)(Y) Xk — vtiY) € I(V) (1 < k < n). In particular, we have that 
the identities 

(2.7) ^(u)X k =v k (u)(l<k<n) 

hold in Q[V]. Observe that the minimal polynomial m u (Y) is square-free, since 
the linear form u separates the points of V. Therefore, m u (Y) and dm u /dY(Y) 
are relatively prime. Thus, multiplying modulo m u (Y) the polynomials Vk(Y) 
by the inverse of (dm u /dY)(Y) modulo m u (Y) we obtain polynomials Wk(Y) := 
(dm u /dY)~ 1 Vk(Y) (1 < k < n) of degree at most D — 1 such that 

(2.8) X k = w k {u) 

holds in Q[V] for 1 < fc < n. The polynomials m u , Wi, . ■ ■ , w n G Q[Y] form a 
geometric solution of V. 

Now, suppose that we are given an algorithm over Q(A) for computing the 
minimal polynomial of the linear form U = AiAi + • • • + A n X n . Suppose further 
that we are given a separating linear form u :— u\X\ + • • • + u n X n G Q[Ai, . . . , X n ] 
such that the vector (ui, . . . ,u n ) does not annihilate any denominator in Q[A] of 
any intermediate result of the algorithm ty. In order to compute the polynomials 
vi, . . . , v n of (|2.7|l . we observe that the Taylor expansion of mjj(A, Y) in powers of 
A — u := (Ai — ui, . . . , A„ — u n ) has the following expression: 

mu (A, Y) = m u (Y) + ]T (^{Y)X k - v k {Y)) (A fc - u k ) mod(A - u) 2 . 

k=l 

We shall compute this first-order Taylor expansion by computing the first-order 
Taylor expansion of each intermediate result in the algorithm ^ . In this way, each 
arithmetic operation in Q(A) arising in the algorithm 'J becomes an arithmetic 
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operation between two polynomials of Q[A] of degree at most 1, and is truncated 
up to order (A — u) 2 . Since the first-order Taylor expansion of an addition, multi- 
plication or division of two polynomials of Q[A] of degree at most 1 requires 0{n) 
arithmetic operations in Q, we have that the whole step requires O(nT) arithmetic 
operations in Q, where T is the number of arithmetic operations in Q(A) performed 
by the algorithm ^. 

Finally, the computation of the polynomials Wi, . . . ,w n of Q2.8jl requires the 
inversion of dm u /dY modulo m u (Y) and the modular multiplication Wk(Y) := 
(dm u /dY)~ 1 Vk(Y) for 1 < k < n. These steps can be executed with additional 
0(nM(Z))) arithmetic operations in Q. Summarizing, we have the following result: 

Lemma 2.4. Suppose that we are given: 

(1) an algorithm ^ in Q(A) which computes the minimal polynomial mu £ 
Q[A, Y] of U :— AXi + ■ ■ ■ + A„A"„ with T arithmetic operations in Q(A), 

(2) a separating linear form u := U\X\ + ■ • ■ + u n X n S Q[-Xi, . . . , X n ] such that 
the vector (u±, . . . , u n ) does not annihilate any denominator in Q[A] of any 
intermediate result of the algorithm ^ . 

Then a geometric solution of the variety V can be (deterministically) computed with 
0(n(T + M(_D))) arithmetic operations in Q. 

3. Statement of the problem and outline of the algorithm 

Let Ai,...,A„ be fixed finite subsets of Z™ with € Aj for 1 < i < n 
and let D := MV(Q\, . . . , Q n ) denote the mixed volume of the polytopes Qi := 
conv( Ai ),..., Q n ■— conv(A n ). Assume that D > holds or, equivalently, that 
dim Q^iei Qi) ^ 1^1 f° r every non-empty subset / C {1, . . . , n} (see, for instance, 
[571 Chapter IV, Proposition 2.3]). 

Let /i, . . . , /„ G Q[X] be polynomials defining a sparse system with respect to 
Ai, . . . , A„ and let d\, . . . , d n be their total degrees. Let d := {d\, . . . , d n }. Suppose 
that fi,...,f n define a zero-dimensional variety V in A" . As in the previous 
section, we group equal supports into s < n distinct supports . . . , .A( s ). Write 
A := (A^, . . . , A^ s ') and denote by ke the number of polynomials fi with support 
A {e) iorl<£<s. 

From now on we assume that we are given a sufficiently generic lifting func- 
tion lo := (u>i, . . . , u> s ) and the fine-mixed subdivision of A induced by u). We 
assume further that the function uii : Ay) — > 1 takes only nonnegative values and 
oje(0, . . . , 0) = for every 1 < £ < s. The lifting function to and the corresponding 
fine-mixed subdivision of A can be used in order to define an appropriate deforma- 
tion of the input system, the so-called polyhedral deformation introduced by Hubcr 
and Sturmfels in [251 - O ur purpose here is to use this polyhedral deformation to 
derive a symbolic probabilistic algorithm which computes a geometric solution of 
the input system f\ = 0, . . . , f n = 0. 

Since the polyhedral deformation requires that the coefficients of the input poly- 
nomials satisfy certain generic conditions, we introduce some auxiliary generic poly- 
nomials gi, . . . ,g n with the same supports Ai, . . . , A„ and consider the perturbed 
polynomial system defined by hi := fi + gi for 1 < i < n. The genericity condi- 
tions underlying the choice of g\, . . . ,g n and hi, ■ . ■ , h n are discussed in Section f3.ll 
We observe that if the coefficients of the input polynomials fi, ■■ ■ ,f n satisfy these 
conditions then our method can be directly applied to fi, . . . , f n . 



12 



G. JERONIMO, G. MATERA, P. SOLERNO, AND A. WAISSBEIN 



Otherwise, we first solve the system hi = 0, . . . , h n = and then recover the 
solutions to the input system fi = 0, . . . , /„ = by considering the homotopy 
h + (1 - T)g x = ■•• = /„ + (1 - T)g n = 0. 

3.1. The polyhedral deformation. This section is devoted to introducing the 
polyhedral deformation of Huber and Sturmfels. 

Let hi :— X) 9 gA c i,qX q for 1 < i < n be polynomials in Q[X] and let V\ denote 
the set of their common zeros in A™. For i = 1, . . . , n, let li be the (unique) integer 
with Aj = A^*', and let u5j := uj£ i be the lifting function associated to the support 
Aj. In order to simplify notations, the n-tuple u := (u>i, . . . , ui n ) will be denoted 
simply by w = (oj\, . . . ,w„). As before, we denote by := {(q,we(q)) G : 
q G C^} the graph of any subset of A™' for 1 < £ < s, and extend this 
notation correspondingly. For a new indeterminate T, we deform the polynomials 
hi, ■ ■ . , h n into polynomials hi, . . . , h n G Q[X, T] defined in the following way: 

(3.1) hi{X,T) := Ci, 9 X*T*M (1 < i < n). 

<?eA ; 

Let / denote the ideal of Q[X, T] generated by hi,...,h n and let J denote the 
Jacobian determinant of hi,...,h n with respect to the variables Xi, . . . , X n . We 
set 

(3.2) V := V(I : J°°) C A" +1 . 

We shall show that, under a generic choice of the coefficients of hi, . . . , h n , the 
system defined by the polynomials in 1)3.1(1 constitutes a deformation of the input 
system hi = 0, . . . , h n = 0, in the sense that the morphism ir : V — > A 1 defined by 
7r(x, t) := t is a dominant map with 7r _1 (l) = V\ x {1}. We shall further exhibit 
degree estimates on the genericity condition underlying such choice of coefficients. 
These estimates will allow us to obtain suitable polynomials hi, . . . , h n by randomly 
choosing their coefficients in an appropriate finite subset of Z. 

According to [2H1 Section 3], the solutions over an algebraic closure Q(T) of 
Q(T) to the system defined by the polynomials 1)3. lfl are algebraic functions of the 
parameter T which can be represented as Puiseux series of the form 

(3.3) x{T) := (xioT 71 / 7 ^ 1 + h-o t., . . . , x n0 T^/^ + h-o t.), 

where 7 :— (71, . . . ,7„,7 n +i) G Z n+1 is an inner normal with positive last coordi- 
nate 7„+i > of a (lower) facet C = (C^ 1 ', . . . ,C^) of A of type (ki, . . . ,k s ), and 
Xo ■= (xio, ■ ■ ■ , x n o) G (C*)™ is a solution to the polynomial system defined by 

(3.4 7 ) := C ^ X<1 (!<*<«), 

where, as defined before, li is the integer with 1 < li < s and Ai = A^'. We shall 
"lift" each solution £0 to this system to a solution of the form 1)3.3(1 to the system 
defined by 1(3. l|l . This means that, on input xq, we shall compute the Puiseux series 
expansion of the corresponding solution ((3.3(1 truncated up to a suitable order. 
Let 

(3.5) V Qn := {x G (C*)" : h^{x) = 0, . . . , h^(x) = 0}. 

A particular feature of the polynomials | |3.4 7 | ) which makes the associated equation 
system "easy to solve" is that the vector of their supports is (C^') kl x ■ ■ • x (C^) fcs , 
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where (C^, . . . , C«) is a cell of type (ki, . . . , k s ) in a fine-mixed subdivision of 
A. Therefore, for every 1 < £ < s, the set C^' consists of ki + 1 points and hence, 
up to monomial multiplication so that each polynomial has a non-zero constant 
term, the (Laurent) polynomials in j3.4 7 [ ) are linear combinations of n + 1 distinct 
monomials in n variables. 

Denote T C Z™ +1 the set of all primitive integer vectors of the form 7 :— 
(71, . . . ,7„,7„+i) G Z™ +1 with 7n+1 > for which there is a cell C = (C*W , . . . , C?W) 
of type (fci, . . . , k s ) of the subdivision of A induced by u such that C has inner nor- 
mal 7. 

Fix a cell C — {C^\ . . . , C^) of type (fci, . . . , of the subdivision of A induced 
by uj associated with a primitive inner normal 7 G T with positive last coordinate. In 
order to lift the points of the variety Vb. 7 of . 5|> to a solution of the system defined 
by the polynomials in 13. If) , we will work with a family of auxiliary polynomials 
/ii i7 , . . . , ft,„ i7 € Q[X, T] which we define as follows: 

(3.6) h i<1 {X,T):=T- m *h i (T^X 1 ,...,T^X n ,T^) (1 < i < n) 

where nij G Z is the lowest power of T appearing in hi(T 11 X\, . . . , T lrl X ni T 7 ™ +1 ) 
for every 1 < i < n. Note that the polynomials obtained by substituting T = 
into /ii. 7 , . . . , /i„ !7 are precisely those introduced in ( |3.4 7 | |. 

3.2. On the genericity of the initial system. Here we discuss the genericity 
conditions underlying the choice of the polynomials g\ , . . . , g n that enable us to 
apply the polyhedral deformation defined by the lifting form u> to the system hi := 
fi + gi = 0, . . . , h n := /„ + g n = 0. 

The first condition we require is that the set of common zeros of the perturbed 
polynomials hi, ■ . ■ , h n is a zero-dimensional variety with the maximum number of 
points for a sparse system with the given structure. More precisely, we require the 
following condition: 

(HI) The set V\ := {x 6 A™ : h\(x) — 0, . . . ,h n (x) = 0} is a zero-dimensional 
variety with D := MV(Qi, . . . , Q n ) distinct points. 

In addition, we need that the system ( |3.4 7 [ ) giving the initial points to our first 
deformation for every 7 G T has as many roots as possible, namely the mixed 
volume of their support vectors. 

For each cell C = (C^\ . . . , C^) of type [hi,..., k s ) of the induced fine-mixed 
subdivision, set an order on the n+ 1 points appearing in any of the sets after 
a suitable translation so that 6 for every 1 < I < s. Assume that G Z™ is 
the last point according to this order. Denote 7 € Z™ +1 the primitive inner normal 
of C with positive last coordinate. Consider the n x (n + 1) matrix whose ith row is 
the coefficient vector of /i|° 7 in the prescribed monomial order and set M. 1 G Q" xn 
and B~f G Q nxl for the submatrices consisting of the first n columns (coefficients of 
non-constant monomials) and the last column (constant coefficients) respectively. 
Then, the coefficients of gi , . . . , g n are to be chosen so that the following condition 
holds: 

(H2) For every 7 G T, the (n x n)-matrix M. 1 is nonsingular and all the entries 
of (M 1 )~ 1 B 1 are nonzero. 

Our next results assert that the above conditions can be met with good proba- 
bility by randomly choosing the coefficients of gi, . . . ,g n in a certain set S C Z. We 
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observe that our estimate on the size of S is not intended to be accurate, but to show 
that the growth of the size of the integers involved in the subsequent computations 
is not likely to create complexity problems. 

Let {Oj i9 : 1 < i < n, g 6 Aj} be a set of new indeterminates over Q. For 
1 < i < n, write f2j := (fij,g : q G Aj) and let Hi G Q[f2j,X] be the generic 
polynomial 

(3.7) Bi(ni,X) := J2 n hi X9 

with support Aj and Ni :— #Aj coefficients. Let fi := (f2i, . . . , f2„) and let N :— 

Nx H h iV n be the total number of indeterminate coefficients. 

We start the analysis of the required generic conditions with the following quan- 
titative version of Bernstein's result on the genericity of zero-dimensional sparse 
systems (see pfl Theorem B], |2H1 Theorem 6.1]): 

Lemma 3.1. There exists a nonzero polynomial P^ G Q[f2] with degP^ < 

3n 2n+l d 2n-l guch that j Qr any c £ qN with p(0)( c ) ^ Q> the system Hi{d,X) 
( ) H n (c n , X) — has D solutions in A 11 , counting multiplicities. 

Proof. Due to 25, Theorem 6.1] combined with 133], the system H\(ci, X) — 0, . . . , 
H n (c n , X) = has D solutions in A™ counting multiplicities if and only if for every 
facet inner normal fi £ Z n of Q\ + • • ■ + Q n , the sparse resultant Res^f 1 ,...^;; does 
not vanish at c := (ci, . . . , c„). Here A^ 1 denotes the set of points of Aj where the 
linear functional induced by /x attains its minimum for 1 < i < n. 

Therefore, the polynomial P^ :— Jl^ ^ es A , 1 '....,AK e where the product 

ranges over all primitive inner normals /i G Z n to facets of Qi + • • • + Q n , satisfies 
the required condition. 

In order to estimate the degree of p(°\ we observe that for every facet inner 
normal /i G Z" the following upper bound holds: 

n 

deg(Res A » ..., AS ) < MV(At, . . . , Af, . . . , A£) < nd"" 1 , 
i=i 

where <i := max{di, . . . , d n }- On the other hand, it is not difficult to see that 
the number of facets of an n-dimensional integer convex polytope P C M™ which 
has an integer point in its interior is bounded by nlvolRn(P). Now, taking P := 
(n + 1)Q, we obtain an integer polytope with the same number of facets as Q 
having an integer interior point. Then, the number of facets of Q is bounded by 
n\vol Rn {P) = n\vol R n((n + l)Q) = (n + l) n n! wJr«(Q) < (n + \) n {nd) n , since Q 
is included in the n-dimensional simplex of size nd. This proves the upper bound 
for the degree P^ of the statement of the lemma. □ 

The next lemma is concerned with the genericity of a smooth sparse system. 

Lemma 3.2. With the same notations as in Lemma \H.l\ and before, there exists a 
nonzero polynomial P^ G Q[f2] of degree at most 4n 2n+1 G? 2 ™ -1 such that for any 
c G Q> N with P^(c) ^ 0, the system Hi(a,X) = 0, . . .,H n (c n ,X) = has exactly 
D distinct solutions in A™. 
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Proof. Consider the incidence variety associated to (Ai, . . . , A„)-sparse systems, 
namely 

W := {(x, c) G (C*) n x {A Nl x • • • x A N ") : ^c^ q x q = for 1 < i < n}. 

As in |4(J1 Proposition 2.3], it follows that W is a Q-irreducible variety. Let tvq : 
W — > A^ 1 x ■ ■ • x A Nn be the canonical projection, which is a dominant map. 

By 37 Chapter V. Corollary (3.2.1)], there is a nonempty Zariski open set 
W(Ai, . . . , A n ) C A^ 1 x ■ • • x A N " of coefficients c = {c\, . . . , c„) for which the 
polynomials Hiipi, X), . . . , H n (c n , X) have supports Ai, . . . , A„ respectively and 
the set of their common zeros in (C*)" is a non-degenerate complete intersection 
variety. Then, the Jacobian Jjj := det(dHi/dXj)i<ij< n does not vanish at any 
point of 7Tn 1 ( c ) f° r ever Y c G U(Ai, . . . , A„). 

Let Q(f2) =— > Q(W) be the finite field extension induced by the dominant projec- 
tion 7tq . By the preceding paragraph we have that the rational function defined by 
Jh in Q(W) is nonzero. Therefore, its primitive minimal polynomial Mj G Q[0, T] 
is well defined and satisfies the degree estimates 

n n 

deg a Mj < deg W ■ deg J H < f](d t + 1) • £ < 2"d n+1 n 

i=i i=i 

(see E3, |H3)- 

Let pW := P {0 Hdf\ where p(°) is the polynomial given by Lemma 13.11 and 
Mj -* denotes the constant term of the expansion of Mj in powers of T. We claim 
that satisfies the requirements of the statement of the lemma. Indeed, let 
ceQ N satisfy pW(c) ^ 0. Then P(°)(c) 7^ holds and so, LemmalO implies that 
H\{c, X) = ■ ■ ■ = H n (c,X) = is a zero-dimensional system. Furthermore, Mj°\c) 
is a nonzero multiple of the product IIa;g 7 r _1 (c) Jh(c,x). Thus, the non-vanishing 

of Mj°\c) shows that all the points of 7Tf 2 1 (c) are smooth and therefore, from e.g. 
|371 IV, Theorem 2.2], it follows that 7r^ 1 (c) consists of exactly D simple points 
in (C*) n . Moreover, combining the assumption that £ A; for 1 < i < n with 
|331 Theorem 2.4], we deduce that 7Tq 1 (c) consists of D simple points in A". The 
estimate degMj 0) < deg Mj < 2 n d n+1 n < ^(n+i^n-i j mp ii cs t k e statement of 
the lemma. □ 

Finally, we exhibit a generic condition on the coefficients hi, . . . ,h n which implies 
that assumption (H2) holds. 

Lemma 3.3. With the previous assumptions and notations, there exists a nonzero 
polynomial P^ 2 ' G Q[^2] with degP*- 2 -* < n(n + l)#r such that for every c :— 
(ci, ...,Cn) G Q N with P {2 \c) ^ 0, the polynomials hi := H^c^X) (1 < i < n) 
satisfy condition (H2). 

Proof. Fix a primitive integer inner normal 7 G V to a lower facet of A. Let 
M. 1 G Q[il] nx ™ and £> 7 G Q[f2]" xl be the matrices constructed from the generic 
polynomials Hi, ... , H n G Q[fi][X] as explained in the paragraph preceding condi- 
tion (H2). Let -Do, 7 G be the (non-zero) determinant of M. 1 , and for every 
1 < i < n i let -Dj',7 be the determinant of the matrix obtained from M. 1 by replac- 
ing its jth column with B 1 . Set P 7 := 11^=0 Af,7- Finally, take P' 2 ' :— n 7 er 
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By Cramer's rule, whenever P^ 2 \c) ^ 0, we have that the system hi, . . . ,h n with 
coefficient vector c = (ci, . . . , c n ) meets condition (H2). 

The degree estimate for follows from the fact that degP 7 < n(n + 1) holds 
for every 7 G T, since each of the entries of the matrices whose determinants are 
involved has degree 1 in the variables f2. □ 

Now, we are ready to state a generic condition on the coefficients of hi, . . . , h n 
which implies that (HI) and (H2) hold. 

Proposition 3.4. Under the previous assumptions and notations, there exists a 
nonzero polynomial P G Q[fi] with degP < 4n 2n+1 <i 2n_1 +n(n+ 1)D such that for 
every c G Q N with P(c) ^ 0, the polynomials hi := Hi(ci,X) (1 < i < n) satisfy 
conditions (HI) and (H2). 

Proof. Set P := p( 1 )p( 2 ) ; where PW is the polynomial of the statement of Lemma 
E3and P( 2 ) is the one defined in the statement of Lemma f3. 31 The result follows 
from Lemmas 13.21 and 13.31 and the upper bound #r < D for the cardinality of 
the set of the distinct inner normal vectors considered (one for each cell of type 
(ki, ■ ■ ■ , k s ) in the given fine-mixed subdivision). □ 

3.3. Outline of the algorithm. Now we have all the tools necessary to give 
an outline of our algorithm for the computation of a geometric solution of the 
(sufficiently generic) sparse system h\ = ■ ■ ■ = h n = 0. 

With notations as in the previous subsections, we assume that a fine-mixed 
subdivision of A induced by a lifting function to is given. This means that we are 
given the set F of inner normals of the lower facets of the convex hull of A, together 
with the corresponding cells of the convex hull of A. In addition, we suppose that 
our input polynomials h\,...,h n G Q[X] satisfy conditions (HI) and (H2) and 
denote by V\ C A" the affine variety defined by hi, . . . , h n . 

First, we choose a generic linear form u G Q[X] such that: 

• u separates the points of the zero-dimensional varieties Vi and Vb,-y for 
every 7 G T. This condition is represented by the nonvanishing of a certain 
nonconstant polynomial of degree at most 2D 2 . 

• An algorithm for the computation of the minimal polynomial of u in Vo, 7 
to be described below can be extended to a computation of a geometric 
solution of Vb,7 according to Lemma T2. 41 for every 7 £ T. This condition is 
represented by the nonvanishing of a nonconstant polynomial of degree at 
most 4D^ for each 7 G T. 

• An algorithm for the computation of the minimal polynomial of u in V to be 
described below can be extended to a computation of a geometric solution 
of V according to Lemma l2~4l This application of Lemma l2~4l rea uires that 
the coefficient vector of the linear form u does not annihilate a nonconstant 
polynomial of degree at most 4D 4 . 

Fix p > 2. From Theorem 12.21 it follows that a linear form u satisfying these 
conditions can be obtained by randomly choosing its coefficients from the set 
{1, . . . , 6/9-D 4 } with error probability at most 1/p. 

Next we compute the monic minimal polynomial fh u G Q(T)[Y] of the linear 
form u in the curve V introduced in (|3.2J) . For this purpose, we approximate the 
Puiseux series expansions of the branches of V lying above by means of a symbolic 
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(Newton-Hensel) "lifting" of the common zeros of the zero-dimensional varieties 
Vb,7 C A" defined by the polynomials 1 |3.4 7 | ) for all 7 6 T (see Section 0J. 

This in turn requires the computation of a geometric solution of Vo. 7 for every 

7 e T. By means of a change of variables we put the system hf^ = ■ ■ • = hn)y = 
defining the variety Vo, 7 into a "diagonal" form (see Subsection 14. II below) . which 
allows us to compute the minimal polynomial mi° 7 of u in Vo, 7 . Since the linear 
form u satisfies condition 2 of the statement of Lemma 12.41 from this procedure we 
derive an algorithm computing a geometric solution of Vo, 7 according to Lemma 
l2~4l 

Then we "lift" this geometric solution to a suitable (non-archimedean) approx- 
imation rhj of a factor m 7 (over Q(T)) of the desired minimal polynomial fh u of 
u. In the next step we obtain the minimal polynomial fh u — H 7 er m i f ronl the 
approximate factors m 7 , namely, we compute the dense representation of the co- 
efficients (in Q(T)) of fh Ul using Pade approximation (see Subsection 14.21 below). 
Finally, we apply the proof of Lemma 12.41 to derive an algorithm for computing a 
geometric solution of the variety V. 

In the last step we compute a geometric solution of the variety V\ by substituting 
1 for T in the polynomials that form the geometric solution of V. 

The whole algorithm for solving the system h\ = • • • = h n = may be briefly 
sketched as follows: 

Algorithm 3.5. 

• Choose the coefficients of a linear form u <E Q[X] at random from the set 
{1,...,6 P # 4 }- 

• For each 7 £ T : 

— Find a geometric solution of the variety Vo l7 defined in l|3.5|l . 

— Obtain a straight-line program for the polynomials /ii. 7 , . . . , h na de- 
fined in 13. 6|) from the coefficients ofhx,...,h n and the entries of 

7 e Z" +1 . 

— "Lift" the computed geometric solution of Vo. 7 to an approximation 
rh 1 of the factor m 7 of m u by means of a symbolic Newton-Hensel 
procedure. 

• Obtain a geometric solution of the curve V : 

— Compute the approximation fh u := Il 7 er ™"i °f ™u- 

— Compute the dense representation of m u from fn u using Pade approx- 
imation. 

— Find a geometric solution ofV applying the proof of Lemma \2.4\ 

• Substitute 1 for T in the polynomials which form the geometric solution 
of V computed in the previous step to obtain a geometric solution of the 
variety V\ . 

4. Solution of the variety V 

4.1. Geometric solutions of the starting varieties. In this subsection we ex- 
hibit an algorithm that computes, for a given inner normal 7 £ T, a geometric 
solution of the variety Vb. 7 C (C*) n defined by the polynomials hf^ (1 < i < n) 
for polynomials h\, . . . , h n satisfying assumptions (HI) and (H2). This algorithm 
is based on the procedure presented in |2"5] . 
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Fix a cell C — [C^\ . . . , C^ 8 ') of type (fci, . . . , k s ) of the given fine-mixed subdi- 
vision of A and let 7 G r be its associated inner normal. For 1 < I < s, we denote 



by hf \ . . . , h)£' the polynomials in the set {h^L . . . , hn,'-y} that are supported in 
C^. In the sequel, whenever there is no risk of confusion we will not write the 
subscript 7 indicating which cell we are considering. 



,(°) 



Our hypotheses imply that h\ 



, h { £ are 



-linear combinations of precisely 



ki + l monomials in Q[X] and, up to a multiplication by a monomial, we may assume 
one of them to be the constant term. Denote these monomials by X ai >° , . . . , X ae - k e , 
with a L0 := £ Z n . Let be the matrix of Qk f x(k e +i) for which thc f u owmg 

equality holds in Q[X, X' 



-life/ 



(4.1) 



M 



if) 



/X aik i\ 




{X a *-° j 





and let MS^ denote the square {kt x fc^)-matrix obtained by deleting the last 
column from M. w . Set 

(M™ •■• \ 
M.W ■■■ 



M 







where here represents different block matrices with all its entries equal to € Q. 
Then M. is the matrix defined by the coefficients of the nonconstant terms of the 
(Laurent) polynomials h^, . . . , hn)y, up to a translation. 

Due to condition (H2) we have that the matrix M. is invertible, which in turn 
implies that the square matrices M.^ are invertible for 1 < I < s. Following 25 , 
we apply Gaussian elimination to the matrix for 1 < ( < s and obtain a set 

of kp + 1 binomials 



(\ 
1 

\o 



ecu* \ ( xae ' ke \ ( Xat '" 1 - c ° 



1 



J \ X"^ 1 J 



\ 



J 



that generate the same linear subspace of Q[X,X 1 ] as the polynomials in l|4.1J) . 
Therefore, for 1 < I < s the set of common zeros in (C*) n of the polynomials 

X "- 1 = c„. ,. Putting 



hi, . • • , hi} is given by the system X a '- k 



these s systems together, we obtain a binomial system of the form 



(4.2) 



X' 



Pi, 



,X a " =p n , 



with a.i G Z n and p, gQ \ {0} (1 < i < n), that defines the variety Vo l7 . Note that 
the second part of condition (H2) ensures the non- vanishing of the constants pi for 
1 < i < n. 

Now, let £ denote the [n x n)-matrix whose columns are the exponent vectors 
ai,...,a„. Using 51, Proposition 8.10], we obtain unimodular matrices K = 

(ki,j)l<i.j<n, l<z ,j <n 

of Z nxn , and a diagonal matrix diag(ri, . . . ,r n ) G 
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Z" xn which give the Smith Normal Form for £, i.e., matrices such that the identity 

(4.3) K-£-L = dmg(n,...,r n ) 
holds in Z" xn . We observe that the upper bound 

(4.4) log||A'|| < (4n + 5)(logn + log||£||) 

holds, where ||A|| denotes the maximum of the absolute value of the entries of a 
given matrix A |511 Proposition 8.10]. 

Let Z\, . . . , Z n be new indeterminates, and write Z := (Z\, . . . , Z n ). We intro- 
duce the change of coordinates given by Xi := z\ 1,% ■ ■ ■ Zn n ' z for 1 < i < n. Making 
this change of coordinates in 14.2fl we obtain the system 

Zi — pi, . . . , /j — Pn, 

which is equivalent to the "diagonal" system 

n n 

Z? = \[{Z Ka J = H/'; =-Qj (1 < 3 < n). 

i=l i=l 

Inverting some of the coefficients qj if necessary we may assume without loss of 
generality that the integers r±, . . . , r n are positive. 

We hrst describe an algorithm for computing a geometric solution of the variety 
Vo. 7 C A™ in the coordinate system of A™ defined by Z\, . . . , Z n . This algorithm 
takes as input the set of polynomials Z^ 1 —qi,..., Z^ n —q n G Q[Zi, . . . , Z n ] defining 
Vo. 7 in the coordinates Z\, . . . , Z ni and outputs a linear form u G Q[Zi, . . . , Z n ] 
which separates the points of Vb, 7 , the minimal polynomial G Q[Y] of u in Vo,7 
and the parametrizations of Z\, . . - ,Z n by the zeros of m^. 

For this purpose, assume that we are given a linear form u := uiZi + - ■ - + u n Z n G 
Q[Zi, . . . , Z n ] which separates the points of Vo i7 . Observe that the fact that u is a 
separating linear form for Vo l7 implies that u.i ^ holds for i = 1, . . . , n. Let Y, Y 
be new indeterminates and let mi, . . . , m„ G Q[Y] be the sequence of polynomials 
defined recursively by: 

(4.5) mi := u^ ri Y ri ~ qil m, := ReSy(u~ r ' {Y -Yp - qi , m^Y)) for 2 < i < n. 

We claim that the polynomial m„ equals (up to scaling by a nonzero element of 
Q) the minimal polynomial m^ G Q[Y] of the coordinate function induced by 
u in the Q-algebra extension Q Q[Vo, 7 ]. Indeed, for every 2 < i < n, the 
polynomial rrii(Y) is a linear combination of u^ ri (Y — Y) ri — qi and mj_i(y) 
over Q[F,y]. Let u {i) := u x Z\ + • • • + UiZi for 1 < i < n. Then, the identity 
u~ Ti (u^ — u^^Y 1 — qi — holds in Q[Vb, 7 ]. Thus, assuming inductively that 
mi-i^" 1 )) = in Q[Vb, 7 ], it follows that mi (u^) = in Q[Vb, 7 ] as well. Taking 
into account that deg m„ < r\ . . . r n and that is a nonzero polynomial of degree 
D-y := n ■ ■ • r n = #(Vo, 7 ), we conclude that our claim holds. 

In order to compute the polynomial mj, we compute the resultants in (|4.5I) . Since 
the resultant ReSy (u~ Ti (Y — Y) Ti — qi, m,_i (Y)) is a polynomial of Q[Y] of degree 
r\ ■ ■ ■ r.i, by univariate interpolation in the variable Y we reduce its computation to 
the computation of ri • • -Ti + 1 resultants of univariate polynomials in Q[y]. This 
interpolation step requires 0(M(rf • ■ • rf )) arithmetic operations in Q and does not 
require any division by a nonconstant polynomial in the coefficients u±, . . . , u n (see, 
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e -g-j |H]> jlOp. Each univariate resultant can be computed using the algorithms in 
e.g. with M(ri • • • r<) arithmetic operations in Q. Altogether, we obtain an 

algorithm for computing the minimal polynomial m% which performs 0(M(Z?^)) 
arithmetic operations in Q. 

Next, we extend this algorithm to an algorithm for computing a geometric solu- 
tion of Vb. 7 as explained in Subsection 12.31 We obtain the following result: 

Proposition 4.1. Suppose that the coefficients of the linear form u are randomly 
chosen in the set {1, . . . , AnpD^}, where p is a fixed positive integer. Then the 
algorithm described above computes a geometric solution of the variety Vq, 7 ( in the 
coordinate system Z\, . . . , Z n ) with error probability at most 1/p using 0(nM(Z3^)) 
arithmetic operations in Q. 

Proof. As proved by our previous arguments, it is clear that the algorithm de- 
scribed computes a geometric solution of Vq~ with the stated number of arithmetic 
operations in Q. There remains to analyze its error probability. 

The only probabilistic step of the algorithm is the choice of the coefficients of 
the linear form u, which must satisfy two requirements. First, u must separate the 
points of the variety Vb, 7 . Since Vb l7 consists of D~ t distinct points of A™, from 
Theorem 12.21 it follows that for a random choice of the coefficients of u in the set 
{1, . . . , 4npZ? 7 }, the linear form u separates the points of Vo, 7 with error probability 
at most l/4np£> 7 < l/2p. 

The second requirement concerns the computation of the univariate resultants 
of the generic versions of the polynomials in 1)4. 5[l . This is required in order to 
extend the algorithm for computing the minimal polynomial mj to an algorithm 
for computing a geometric solution of the variety Vb l7 . We use a fast algorithm 
for computing resultants over Q(A) based on the Extended Euclidean Algorithm 
(EEA for short). We shall perform the EEA over the ring of power series Q[A — u\, 
truncating all the intermediate results up to order 2. Therefore, the choice of 
the coefficients of u must guarantee that all the elements of Q[A] which have to 
be inverted during the execution of the EEA are invertible elements of the ring 
Q[A-SJ. 

For this purpose, we observe that, similarly to the proof of |56l Theorem 6.52], 
one deduces that all the denominators of the elements of Q(A) arising during the 
application of the EEA to the generic version of the polynomials u~ Ti (a—u( l ~ 1 ') ri — 
qi and TOi_i(u^~ 1 * ) ) are divisors of at most r± ■ ■ ■ rj_i polynomials of Q[A] of degree 
2ri ■ • • for any a G Q. This EEA step must be executed for n • • ■ ri distinct 
values of a G Q, in order to perform the interpolation step. Hence the product 
of the denominators arising during all the applications of the EEA has degree at 
most 2nD^. Therefore, from Theorem 12.21 we conclude that for a random choice 
of its coefficients in the set {1, . . . ,4np_D^}, the linear form u satisfies our second 
requirement with error probability at most 1/2,0. 

The lemma follows putting both error probability estimates together. □ 

Finally, we compute a geometric solution of the variety Vo, 7 in the original co- 
ordinate system defined by X\, . . . , X n . 

For this purpose, we compute the minimal polynomial m u G Q[Y] of a linear form 
u = uiXi + • • • + U n X n G Q[Xi, . . .,X n ] in Vb, 7 . Let V , 7 := {x^ 7) , . . . ,x^' 7) }. 
Then we have m u (Y) = Y\f=i{Y — u(xq J ' 7 ')). In order to compute m u , we use 
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the polynomials mj, wi, . . . , w n which form the previously computed geometric 

k (~r) fc (T) 

solution of Vo i7 in the variables Z\, . . . , Z n : from the identities Xj := Z^' % ■ ■ ■ Z^' % 
(1 < i < n) we deduce that m u equals the minimal polynomial of the image of 

the projection ?7„ : Vo, 7 — » A 1 defined by r)^\zi, . . . , z„) := Yl7=i u i z i 1A ' ' ' z n n,i ■ 
Now, the identities Zi = z«i(u), which hold in Q[Vb. 7 ] for 1 < i < n, imply that 

(4.6) n = > ,Uj(wi(n)) M ■■■(w n (u)) "' l 

i=l 

holds in Q[Vb, 7 ], from which we easily conclude that m u satisfies the following 
identity: 

(4.7) m u [Y) - ReSy{Y - ^u^w^Y )) k ^ ■ ■ ■ {w n {Y)) k ^ , m 2 (F)) . 

i=l 

We compute the monomials (wi(u)) 11 • ■ • (w n (u)) ™. ; (1 < i < n) in the right- 
hand side of (|4.6|) modulo mjfK), with 0(n 2 log(max, i:) - |A;^|)M(£) 7 )) additional 
arithmetic operations in Q. From 14.411 it follows that 

0(n 2 log(max \kM\)M{D 7 )) - 0(n 3 log(n||£ 7 ||)M(D 7 )), 

where f 7 is the matrix of the exponents of the cell corresponding to the inner 
normal 7. Observe that all these steps are independent of the coefficients of the 
linear form u we are considering and therefore do not introduce any division by a 
nonconstant polynomial in the coefficients u\, . . . ,u n . 

In the next step we compute the right-hand side of l|4.6|l modulo rriu(Y), with 
0(n_D 7 ) arithmetic operations in Q. Then we compute the resultant l|4.7|l by 
a process which interpolates (|4.7|l in the variable Y to reduce the question to the 
computation of D 7 + l univariate resultants, in the same way as for the computation 
of the resultants in (|4.5I) . This requires 0(M(Z? 7 ) 2 ) arithmetic operations in Q. 

If the linear form u separates the points of Vb )7 , then we can extend the algorithm 
for computing m u (Y) to an algorithm for computing a geometric solution of Vo )7 
with the algorithm underlying the proof of Lemma l2.4l This extension requires that 
the coefficients u±, . . . , u n of the linear form u do not annihilate the denominators 
in Q[A] which arise from the application of the algorithm described above to the 
generic version K\X\ + • • ■ + A n X n of the linear form u. Such denominators arise 
only during the computation of the generic version of the resultant (|4.7I) . Hence, 
with a similar analysis as in the proof of Proposition 14.11 we conclude that, if 
the coefficients of u are chosen randomly in the set {1, . . . ,4pD^}, then the error 
probability of our algorithm is bounded by In conclusion, we have: 

Proposition 4.2. Suppose that we are given a geometric solution of Vb. 7 in the 
coordinate system Z±,...,Z n , as provided by the algorithm underlying Proposi- 
tion \4-l\ and the coefficients of the linear form u are randomly chosen in the set 
{1, . . . , 4pD^}, where p is a fixed positive integer. Then the algorithm described 
above computes a geometric solution of the variety Vo. 7 with error probability at 
most 1/ p using 0(n? log(n||£ 7 ||)M(_D 7 ) 2 ) arithmetic operations in Q. 

Finally, from Propositions 14.11 and 14.21 and the fact that ||-E 7 || is bounded by 
Q := 2 maxi<i< ra {||(7||; q £ A,}, we immediately deduce the following result: 
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Theorem 4.3. Suppose that the coefficients of the linear forms u and u of the state- 
ment of Propositions \4-l] and \4-S\ are chosen at random in the set {1, . . . , 4npD 3 }, 
where p is a fixed positive integer. Then the algorithm underlying Propositions \4-l\ 
and \4-S\ computes a geometric solution of the varieties Vo l7 for all 7 £ T with error 
probability at most 2/p using 0(n 3 log(nQ)M(£>) 2 ) arithmetic operations in Q. 

4.2. The computation of a geometric solution of the first deformation. 

The second step of our algorithm consists in the computation of a geometric solution 
of the curve V of (|3.2J) . This will be done by "lifting" the geometric solutions of 
the varieties Vo i7 computed in the previous section for all 7 G T. 

We recall the definition of the variety V. Let / denote the ideal of Q[X, T] gener- 
ated by the polynomials hi, . . . , h n of which form the polyhedral deformation 
of the generic polynomials hi, . . . ,h n , and let J denote the Jacobian determinant of 
hi, . . . , h n with respect to the variables Xi, . . . , X n . Let V(I) be the set of common 
zeros in A" +1 of hi, . . . ,h n . Then V := V(I : J°°). 

Alternatively, let 7r : V(I) — » A 1 be the linear projection defined by it{x, t) — t. 
Consider the decomposition of V(I) into its irreducible components V{I) = Ul=i 
Suppose that the restriction 7r|c ; : Ci — > A 1 of the projection ir is dominant for 
1 < i < r and is not dominant for r + 1 < i < s. We shall show that V :— Ul=i 
holds, i.e., V is the union of all the irreducible components of V(I) which project 
dominantly over A 1 . Furthermore, we shall show that V C A n+1 is a curve which 
constitutes a suitable deformation of the variety defined by the system hi = ■ ■ ■ = 
h n = 0. For this purpose, we shall use the following technical lemma: 

Lemma 4.4. Let F x ,...,F n £ Q[X,T] and V := {{x,t) E A n+1 : = 
0, . . . , F n (x, t) = 0}. Set I := {F±, . . . , F n ) C Q[X, T] and let J denote the Jacobian 
determinant of Fi , . . . , F n with respect to the variables X. Consider the linear 
projection ir : V — ► A 1 defined by n(x, t) :— t. Assume that #7i-- 1 (<) < D holds 
for generic values of t € A 1 and that there exists a point to G A 1 such that the 
fiber 7r _1 (io) is a zero- dimensional variety of degree D with J(x,to) ^ for every 

(x,t ) e tt" 1 ^,)). 

Let Vdom be the union of all the irreducible components CofV with ir(C) — A 1 . 
Then: 

• Vdom is a nonempty equidimensional variety of dimension 1 . 

• Vdom is the union of all the irreducible components ofV having a non-empty 
intersection with 7r -1 (to)- 

• V dom = V{I: J°°). 

• The restriction 7r|y dom : Vdom — > A 1 is a dominant map of degree D. 

Proof. First we observe that dim(C) > 1 for each irreducible component C of V, 
since V is defined by n polynomials in an (n + l)-dimensional space. 

Let C be an irreducible component of V for which 7r -1 (io)nC 7^ holds. Consider 
the restriction n\c : C — > A 1 of the projection map tt. Then we have that tt\q 1 (to) is 
a nonempty zero-dimensional variety, which implies that the generic fiber of n\c is 
either zero-dimensional or empty. Since dim(C) > 1, the Theorem on the Dimension 
of Fibers implies that dim(C) = 1 and that ir\c ■ C — * A 1 is a dominant map with 
generically-finite fibers. This shows that C C Vdom and, in particular, that Vdom is 
nonempty. 
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Conversely, we have that 7r _1 (io) H C ^ holds for any irreducible component C 
of Vdom- Indeed, assume on the contrary the existence of an irreducible component 
Co not satisfying this condition. Then, there is a point t\ G A 1 having a finite 
fiber 7r" 1 (^i) such that 7r| t 7 1 (ti) and ^Ic^ 1 (^i) have maximal cardinality for every C 
with C n 7T _1 (t ) ^ 0- This implies that #7r" 1 (ti) > #7r -1 (io) = A leading to a 
contradiction. 

We conclude that Vdom is the nonempty equidimensional variety of dimension 
1 which consists of all the irreducible components C of V with 7r _1 (io) n C 7^ 0. 
Furthermore, this shows that the restriction 7r|y dom : Vdom —> A 1 is a dominant 
map of degree D. 

Finally we show that the identity Vd om = V(I : J°°) holds. First, note that the 
irreducible components of V(I : J°°) are all the irreducible components of V where 
the Jacobian J does not vanish identically. Thus, it is clear that Vdom C V(J : J°°), 
since J does not vanish at the points of 7r _1 (io) nC for each irreducible component 
C of Vdom- On the other hand, if C is an irreducible component of V for which 
the projection tt| c : C -> A 1 is not dominant, then C is the set of common zeros 
of the polynomials Fx, . . . , F n ,T — tc for some value tc- Since dim(C) > 1, we 
have that the Jacobian matrix d(Fi, . . . , F n , T — tc)/d(Xx, ■ ■ . , X n , T) is singular 
at every point (x,tc) of C. Hence, its determinant, which equals J, vanishes over 
C. □ 



Now we return to the study of the variety V and show that the assumptions of 
Lemma [4.41 hold. Observe that tt^ 1 ^) = Vt x {t} holds for every ( G A 1 , where 
V t := {x e A™ : hi(x,t) — 0, . . . , h n (x, t) = 0}. Furthermore, the polynomi- 
als hi(X,t), . . . ,h n (X,t) are obtained by a suitable substitution of the variables 
57 of the generic polynomials Hi, ... , H n e X] with supports Ai, . . . , A n 

introduced in l|3.7|l . Indeed, if c = (ci,...,c„) is the vector of coefficients of 
h\, . . . ,h n , the coefficient vector of hi(X,t) (1 < i < n) is (cj g^^^gA* for every 
( 6 A 1 . By Lemma f3.ll there exists a nonzero polynomial p(°) G such that 

for any c' = (c[, . . . ,c' n ) with p(°^(c') 7^ the associated sparse system defines 
a zero-dimensional variety. In particular, the coefficients c = (ci,...,c„) of our 
input polynomials hi := Hi(ci,X), . . . , h n = H n (c n ,X) satisfy P(°'(c) ^= 0. This 
shows that the polynomial G Q[T] obtained by substituting Oj i9 1— » Cj g T w< ^^ 
(1 < i < n, q G Aj) in the polynomial P^ ^ is nonzero, since it does not vanish at 
T = 1. We conclude that V* is a zero-dimensional variety for all but a finite number 
of t G A 1 . Thus, 7r _1 (t) is finite for generic values of t G A 1 . 

Finally, by condition (HI), the fiber 7r _1 (l) = V(hi, . . . , h n ) x {1} is a zero- 
dimensional variety of degree D = deg(7r) and the Jacobian determinant J := 
dct(dhi/dXj)i<ij< n does not vanish at any of its points. On the other hand, the 
fact that #7r -1 (i) < D holds for generic values t G A 1 follows from the BKK 
theorem. 

This shows that the variety V(I) and its defining polynomials hi, . . . , h n satisfy 
all the assumptions of Lemma 14.41 Thus, we have: 

Lemma 4.5. The variety V C A n+1 is a curve. Furthermore, every irreducible 
component ofV has a nonempty intersection with the fiber 7r _1 (l) of the projection 
map 7r : V — > A 1 . 
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4.2.1. Generic linear projections of V. In order to compute a geometric solution 
of the space curve V, we shall first exhibit a procedure for computing the minimal 
polynomial of a generic linear projection of V. Let u 6 Q[Xi, . . . , X n ] be a linear 
form which separates the points of the "initial varieties" Vo. 7 for all the inner 
normals 7 :— (71, ... , 7n+i) of the lower facets of the polyhedral deformation under 
consideration. Let ir u : V — > A 2 be the morphism defined by ir u (x,t) :— (t,u(x)). 
Since the projection map ir : V — > A 1 defined by n(x, t) := t is dominant, it follows 
that the Zariski closure of the image of ir u is a Q-definable hypersurface of A 2 . 
Denote by M u € Q[T, Y] a minimal defining polynomial for this hypersurface. For 
the sake of the argument, we shall assume further that the identity deg(7r) = D, 
and thus deg y M u = D, hold. 

We can apply estimate l|2.2|) of Lemma 12.31 in order to estimate deg T M u in 
combinatorial terms (compare with (4HI Theorem 1.1]). Indeed, let Qi, . . . ,Q n C 
R" +1 be the Newton polytopes of the polynomials h\, . . . , h n of l|3.1|) . and let A C 
R" +1 be the standard unitary simplex in the plane {T = 0}. Then the following 
estimate holds: 

(4.8) deg T M u <E:= MV n+1 (A, Q u . . . , Q n ). 

Furthermore, equality holds in l|4.8|l for a generic choice of the coefficients of the 
polynomials hi and the linear form u. 

More precisely, we shall exhibit a procedure for computing the unique monic 
multiple in Q(T)[F] of M u of degree D. This polynomial can be alternatively 
defined as explained in what follows: 

Since the projection map 7r : V — > A 1 is dominant, it induces an extension 
Q[T] <^-> Q[V], where Q[V] denotes the coordinate ring of V. This variety being 
a curve, Q[V] turns out to be a finitely generated Q[T]-module. Thus, tensoring 
with Q(T), we deduce that Q[V] ®Q(T) is a Q(T)-vector space of finite dimension. 
We claim that Q[V] ® Q(T) = Q[V(I)} ® Q(T) holds. Indeed, since V consists of 
the irreducible components of V(J) which are mapped dominantly onto A 1 by the 
projection 7T, for each of the remaining irreducible components C of V(I), the set 
7r(C) C C is a zero-dimensional Q-definable variety. This implies that 1(C) f)Q[T] ^ 
{0} holds. 

Let rh u be the minimal polynomial of u in the extension Q(T) ^ Q[V] ® Q(T). 
The fact that Q[V] ®Q(T) is finite-dimensional Q(T)-vector space shows that the 
affine variety V :— {x e A"(Q(T)*) : = 0, . . . , /i n (S) = 0} has dimension zero. 

Here Q(T)* := U gGN denotes the field of Puiseux series in the variable 

T over Q (see, e.g., |S7]) and h\,...,h n are considered as elements of Q(T)[X], 
Our hypotheses imply that there exist D distinct n-tuples := (x[ , . . . , Xn ^) € 
(Q(T)*) n of Puiseux series such that the following equalities hold in Q(T)* for 
1 < £ < D: 

(4.9) fti(^,T) = 0, ... , ft n (a;W s T) = 

(see |23). Since Q[V] ® Q(T) is the coordinate ring of the (Q(T)-variety ¥, from 
H4.9|) we deduce that the dimension of Q[V]®Q(T) over Q(T) equals D . Moreover, 
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since as a consequence of our assumptions degy fh u = D holds, we conclude that 

D 

(4.10) fh u = l[(Y-u(x^)). 

i=\ 

Since M U (T, u(X)) 6 I(V), it follows that M U (T, u(X)) = holds in Q[V] ® Q(T), 
from which we conclude that M u is a multiple of fh u by a factor in Q(T)[Y}. Taking 
into account that both are polynomials of degree D in the variable Y and that fh u 
is monic in this variable, we deduce that fh u is the quotient of M u by its leading 
coefficient. 

4.2.2. A procedure for computing fh u . Now we exhibit a procedure for computing 
the minimal polynomial m u , which is based on the expression (|4.10() of rh u in terms 
of the Puiseux expansions (|3 .31) . Then we will apply Lemma l2.4l to this procedure 
in order to obtain an algorithm for computing a geometric solution of the curve V . 

With notations as in Section Em let V C Z" +1 be the set of primitive integer 
vectors of the form 7 := (7 1; . . . ,7 n ,7„+i) G Z n+1 with 7„ + i > for which there 
is a cell C = (C^\ . . . , C^) of type [ki, . . . , k s ) of the subdivision of A induced 
by u> such that C has inner normal 7. As asserted in Section mi if 7 6 T is the 
inner normal of the lifting C of a cell C of type (hi,..., k s ), there exist D 1 := 
fci! • • • k s \ ■ vol(C) vectors of Puiseux series := {x[ 3 ' 7 \ x^) G A n (Q(T)*) 

(1 <i < D y ) of the form 

T 0'.7) ._ \" r tt?)T-Sr 

m>0 

satisfying (|4.9|l . Considering the projection of the branches of V parametrized by 
the -D 7 vectors of Puiseux series x^^ for each 7 £ T, we obtain the following 
element m 7 of Q((T 1/7 " +1 ))[y]: 

£» 7 

(4.11) m 7 := ]J (F-u(a; (j '' 7) )). 

From l|2.1|l we conclude that (|4.1()[) may be expressed in the following way: 

(4.12) m u =Y[ m i- 

Since fh u belongs to Q(T)[Y"] and its primitive multiple M u g Q[T, Y) satisfies 
the degree estimate deg T M u < E, in order to compute the dense representation 
of fh u we shall compute the Puiseux expansions of the coefficients of the factors 
m 7 G Q((T 1 / ln + 1 ])[Y] of fh u truncated up to order 2E. Using Pade approximation 
it is possible to recover the dense representation of fh u from this data. 

Fix 7 G T and set x^' 7 ^ := (a^ m > ■ ■ ■ > ^n'm) for every m > and 1 < j < _D 7 . 
Since 

hi( E sS?^, • ■ ■ , E &^,T) = 

m>0 m>0 
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holds for 1 < j < D-y and 1 < i < n, we have 

(4-13) = T~ m *hi(T^ E m >o 4ff T"*, • • • > T7 " E ra >o a^T"*, T 7 ^ 1 ) 

- /) fV x 0'.7) T m jA 

— '%,7^ Z^m>0 Xm 1 

according to l|3.6|l . Therefore the polynomial m 7 (T 7 ™ +1 ,F) G Q((T))[Y] can be 
expressed in terms of the power series solutions cr^ 7 ' := (<r^ , . ■ ■ , Cn ) '■— 
Em>o x ™ 7)T ™ (! < i < ^7) 01 ft l,7> ■ ■ ■ , ^n, 7 - Indeed, from @JT) it follows that 

m 7 (T 7 ^, Y) = nfi ^ - ELi "i E ro >o *S^ 7i+m ) 



(4.14) 



life - E m >o E? =1 t^T*!™) 
life (^-E m >o"7(x^ 7) )T™) 
nfji(^-%(E m >o^' 7) T" 1 )) =: m» 7 (T,y), 



where u 7 := Ei=i ^i 7 " 17 *-^- We conclude that the Laurent polynomial 
m 7 (T 7ra+1 , Y) G Q((T))[y] may be considered as the minimal polynomial m u (T,Y) 
of the projection induced by u 7 on the subvariety V 1 of A n (Q(T)*) consisting of 
the set of power series {o^ 1,7 - 1 , . . . , er^' 7 ")}. This remark will allow us to compute 
a suitable approximation to the Laurent polynomial m 7 (T 7 ™ +1 , Y) in Q((T))[Y]. 

In order to describe this approximation, we introduce the following terminology: 
for G, G G Q({T)) and s G Z, we say that G approximates G with precision s in 
Q((T)) if the Laurent series G — G has order at least s + 1 in T. We shall use the 
notation G = G mod (T s+1 ). Furthermore, if G, G are two elements of a polynomial 
ring Q((T))[Y], we say that G approximates G with precision s if every coefficient 
a G Q((T)) of G approximates the corresponding coefficient a G Q((T)) of G with 
precision s (in the sense of the previous definition). 

Proposition 4.6. Fix 7 := (71, . . . ,7„) G T and assume that a geometric solution 
of the variety Vo, 7 is given, as provided by Theorem \4-3[ Assume further that the 
coefficients of the linear form u of the given geometric solution ofVo^ are randomly 
chosen in the set {1, . . . , 4pD^} for a given p G N. Then there is an algorithm 
which computes an approximation to the polynomial m Uy G Q((T))[Y] with precision 
2Ej n+1 . The procedure requires 0((nL 7 + n n )M(D 7 )(M(M 7 )M(L> 7 ) + Ej n+1 )) 
arithmetic operations in Q , where M 7 := max{7i, . . . ,7„} and L 7 is the number 
of arithmetic operations required to evaluate the polynomials h% n of and has 

error probability at most 2/ p. 

Proof. Let notations and assumptions be as before. In order to compute the re- 
quired approximation of the polynomial m Ui we first compute the corresponding 
approximation of the polynomials that form a geometric solution of the variety 
K 7 := {cr(J' 7 ); 1 < j < D-y}. Observe that 

{^• 7 >(0);1 < j < D,} = {x^;l<j<D y } 

= v(h^,...,hM)n(c*r 

= v(h ln {x, o), . . . , h n , 7 (x, o)) n (c*)" = v 0n 
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holds. Since det(dh in (X 7 0)/dXk)i<i l k<n(xo' t ) ^ holds for 1 < j < Z? 7 , we may 
apply of the global Newton iterator of |2J (see also [SD] ) in order to "lift" the given 
geometric solution of Vb l7 to the geometric solution of the variety V 1 associated to 
the linear form u £ Q[X] with any prescribed precision. 

Denote ml° 7 , 7 , ■ • ■ , w£)n,~f 6 Q[Y] the polynomials which form the given 
geometric solution of Vb, 7 , as provided by the algorithm underlying Theorem 14.31 
Recall that ml° 7 (u(x[, 3) )) = and (x^' 7) )i = w { °}„ (u(x { Q j) )) holds for 1 < i < n 
and 1 < j < D-y. The global Newton iterator is a recursive procedure whose fcth 

step computes approximations ?7il fe 7 , , . . . , wl!%^ € Q[T, Y] of the polynomials 

™«, 7 , w M ,i j7 , . . . , w u>nn which form the geometric solution of V-y associated with the 
linear form u with precision 2 k for any k > 0. 

Assume without loss of generality that ji > and = min{7i, . . . , 7„} hold for 
1 < i < n. Indeed, if there exists ji < 0, setting 7j := min{7i, . . . , 7„} we have 

D 7 n 

n T ^-° ( T7, ° y - e e ^™ )r7i+m ) 

j — 1 2=1 m>0 

£>7 n 

n(^-^ 7<o E^E^ 7)T7<+m ) 

j — 1 2=1 m>0 

n( y -E^E^ )T7 -" 7iD+m )- 

j — 1 i— 1 m>0 

Since 7^ — 7i > holds for 1 < i < n, this shows that the computation of an 
approximation m u := m 7 (T 7 ™ +1 , Y) can be easily reduced to a situation in which 
7j > holds for 1 < i < n. 

Note that the global Newton iterator cannot be directly applied in order to 
compute the geometric solution of {c^' 7 - 1 ; 1 < j < Z3 7 } associated with the linear 
form u 1 € Q[T][X], because the coefficients of u 7 are nonconstant polynomials of 
Q[T]. Indeed, two critical problems arise: 

(1) Although by hypothesis m 7 separates the points of V-y, it might not separate 
the points of Vo. 7 and it is not clear from which precision on, the corre- 
sponding approximations of the points of Vy are separated by u 7 . Requiring 
Uj to be a separating form for all the approximations of the points of ~V 1 is 
an essential hypothesis for the iterator of [2] which cannot be suppressed 
without causing a significant growth of the complexity of the procedure (see 

EDI, EH). 

(2) The iterator of [21] makes critical use of the fact that the coefficients of 
the linear form under consideration are elements of Q in order to determine 
how a given precision can be achieved. 

Nevertheless, we shall exhibit a modification of the procedure which computes an 
approximation of m u ^(T,Y) with precision 2"f n+ iE without changing the asymp- 
totic number of arithmetic operations performed. 

In order to circumvent Q we require an additional generic condition to be 
satisfied by the coefficients ui, . . . ,u n defining w 7 :— X)"=i UiT li Xi. Recall that 
= E m > (E"=i ^ m 7 l 7 jT™ for every 1 < j < D ? where xg^ := 
for m < 7j. To state this condition, we need the following claim: 
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Claim. Set M 7 := max{'yi, . . . /y n } and let Ai,...,A n be indeterminates over 
C[T,X]. Then the following inequality holds for every 1 < j, h < Dj with j =/= h: 

My n M*y n 

E (X x - r - ■)'"' - E (EM^K- 

m—0 i—1 rn—O i—1 

Proof of Claim. Suppose on the contrary that there exist j ^ h such that 
YIU (£?=i A, ^l 7i )T- = A, z&ljT™ Substituting T^A, 

for A. ( in this identity for i = 1, . . . , n, we have £m=o £™=i Aj x ^m- lz Tm ~ lz = 
YZU Eti A, 4t 7 A ri T m -' w , that is 

E E M^ m = E E a.^t™. 

i—1 m— i—1 m—0 

Substituting for T in this identity, we deduce that 

n n 

VA'i ftT) - V A-a^' 1 ' 7 ' 

/ ^-' L t J 'i.O — / ^H-^-i.o ; 
i=l i=l 

which contradicts the fact that the vectors Xq ' = {x^q\ ■ ■ ■ > x n'o^) (1 — i — ^7) 
are all distinct. This finishes the proof of the claim. 

By the claim we see that the polynomial J2m=o ( Z)?=i A i (^i^m- 7 , _x i,m 7 -7 4 
of Q[A][T] is nonzero, and therefore has a nonzero coefficient ajji G C[A] for every 
1 < j < h < D 1 . Consider the polynomial A 7 (A) := Yli<j<h<D cij,h £ C[A]. Since 
aj t h has degree 1 for every 1 < j < h < £> 7 , it follows that A has degree f^ 7 )- Fur- 
thermore, for every (u\, . . . , u„) 6 C™ with A 7 (mi, . . . , u„) 7^ 0, the corresponding 

polynomial u 7 := £™ = i u i^" 1i ^-i separates the initial terms £m=o x m' 7 ^ m °f the 
power series a* 7 ' (1 < j < Z5 7 ). 

From Theorem 12 .21 we see that for a random choice of the coefficients u\, ■ ■ ■ , u n 
in the set {1, . . . , pD^} the linear form u 7 separates the first M 7 terms of the points 
of V~f with probability at least 1 — 1/p. From now on we assume that it 7 satisfies 
this requirement. 

The algorithm proceeds in three steps. First, it computes a suitable approxima- 
tion to the geometric solution of V 1 associated to the linear form u :— £™ =1 u iXi 
by means of n := |~log(M 7 + 1)] steps of the global Newton iterator of [2"T] . 
This approximation is used in order to obtain the corresponding approximation 
mu 7 °\^ K °i, • • • j Wu^li of the polynomials that form the geometric solution of V~ ( 
associated with u 7 . Finally, we apply an adaptation of the global Newton iterator 
which takes as input the polynomials of the previous step m,u°\ w ^a°\i ■ • • > u4 7 °n 
and outputs the required approximation to the polynomials m u , w u 1, . . . ,w u n 
that form the geometric solution of V 1 associated with u 7 . 

Now we consider the three steps above in detail. The first step takes as in- 
put the given geometric solution ml° 7 , , . . . , u4°n, 7 of Vb, 7 , and performs 
Ko := [~log(M 7 + 1)] times the global Newton iterator of |2J to obtain polyno- 

(«o) 
",1,7 ' 



mials rriu,° , w„ ° ' , . . . , w^n^ £ Q[T, Y] such that the following conditions hold: 



{i)u, Ko de Sy = D ~< and degyml^- 1 < M 7 , 

(«)«,«□ deg F w^- l < D y and deg T w^°' < M 7 for 1 < i < n, 
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{iii)u,. = nfJi (Y - <P«o) mod (T M " I+1 ), 

(iv) u , Ko ap' 7) = (T, ^J 7) ) mod for 1 < i < n. 

Here </?k ' 7 ' ) is the Taylor expansion of order 2 K ° of the power series u(cr^' 7 '), that 
is, ^ 7) := Erio"( x - 7) ) Tm fOT 1 < 3 < D i- 

According to 21, Proposition 7], it follows that this step requires performing 
0((nX 7 + n n )M(D 1 )M(M 1 )) arithmetic operations in Q, where L 7 denotes the 
number of arithmetic operations in Q required to evaluate the polynomials hi^ of 

6(1 . Furthermore, in view of the application of Lemma 12.41 it is important to 
remark that this step does not involve any division by a nonconstant polynomial in 
the coefficients u\, . . . , u n . 

Next we discuss the second step. In this step we obtain approximations 
Tn[^°\ w^\, . . . , w^°n of the polynomials that form the geometric solution of V 1 
associated with u 1 with precision 2 K ° > M~, namely 

• deg y m!^ o) = Dj and deg T m^ < M 7 , 

• deg Y w { *°J < D 1 and deg T w ( *°J < M 7 for 1 < i < n, 
. ml K ; ) = nfja (Y - <^' 7) ) mod (T^+i), 

. <7^' 7) ee ^(T, ^ ,7) ) mod (T M - +1 ) for 1 < i < n. 

Here 0k o is the Taylor expansion of <j>^^ :— « 7 (cr^' 7 )) of order 2 K ° for 1 < j < 
D 1 . 

From conditions (i) UlKo - {iv) UlK0 and the elementary properties of the resultant 
it is easy to see that rrSu^ satisfies the following identity: 

n 

(4.15) m^\Y) = Res Y {Y -J2^^i%(n "»<#(?)). 

The resultant of the right-hand side is computed mod (T M ~< +1 ) by interpolation 
in the variable Y to reduce the problem to the computation of -D 7 resultants, as 
explained in the computation of the resultant in l|4.7|l . These D 1 resultants involve 
two polynomials of Q[T, Y] of degree in Y bounded by £> 7 and are computed mod 
(T m t +1 ). Hence we deduce that this step requires 0(M(D 7 )Z? 7 M(M 7 )) arithmetic 
operations in Q. 

We apply Lemma 12.41 in order to extend this procedure to an algorithm comput- 
ing m^°\ w^°\, . . . ,w^°n. For this purpose, we observe that a similar argument 
as in the proof of Proposition ^. ll proves that the denominators in Q[A] which arise 
during the computation of the -D 7 resultants required to compute the minimal poly- 
nomial of the generic version Y^i=i AjT^-Xi of the linear form u 7 are divisors of a 
polynomial of Q[A] of degree at most 4Z? 7 . Applying Theorem 12. 21 we see that for a 
random choice of the coefficients ui, . . . ,u n in the set {1, . . . , 4pZ? 7 } none of these 
denominators are annihilated with probability at least 1 — 1/p. 

Finally, we consider the third step of the algorithm. For ki := |~log(27„ + ii?+l)] , 
we apply «i — kq times an adaptation of the global Newton iterator of ;2*p| to the 
polynomials m^°\ w^°\, . . . , wC°} n computed in the previous step. In the kth 
iteration step, we compute polynomials mi 7 ,wf£j 1; . . . , w£j, n satisfying: 

• deg Y mi*; = D and deg T m£} < 2 k , 
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• deg Y w^u \ < D and deg T w^\ 1 < 2 k for 1 < i < n, 

• a\ jn) =^(T,^ 7) ) modir 2 ^ 1 ) for 1 < i < n. 

Here <t>£^ is the Taylor expansion of (/fin) := u^{a^'^) of order 2 k for 1 < j < D~ ( . 
In particular, it follows that m Ul is the required approximation to m u with 
precision 2j n+ iE. 

Fix kq < k < K\. We briefly describe how we can obtain an approximation with 
precision 2 k of the polynomials that form the geometric solution of Vy associated 
to the linear form u 7 from an approximation with precision 2 fe ~ 1 . Similarly to 
EH, set A k (T,Y) := u 7 (w$) ~ u^w^) = u 7 (w^) - Y, where tuJJJ is the 
result of applying a "classical Newton step" to w^T^, as described in [21]. Fur- 
thermore, write A m (T, Y) T" 1 " 2 * -1 - m£ -1) ). Since mi^Y + A fe ) = 
mod (T 2k+1 ,m ( u~ 1) ) holds (see [H §4.2]), it follows that 

= mW(y+ A fe ) ee m^Yf A fe ) + T 2 ^ 1+1 A m (r+ A,) mod (T 2fe+1 , m^" 1 )) 

A.^MF) +T 2fc - 1 + 1 A m( r) mod (T* h +\m£-V). 
We conclude that the following congruence relation holds: 

(4.16) <) ^ mfrV - (Afe^^ mod mgp 1 ') mod (T 2 ^ 1 ). 
A similar argument proves the following congruence relation 

dw (k ~ 1} 

(4.17) ^.^^-(A,^ modm^)) mod (T 2 ^ 1 ) for 1 < i < n. 

Each iteration of our adaptation of the global Newton iteration is based on (|4.16(l 
and 1|4.17J) . which are extensions of the corresponding congruence relations of |21j . 

~(k) 

We first compute u>« T by a standard Ncwton-Hcnscl lifting, and then evaluate the 
expressions (|4.16|) and l|4.17|l . With a similar analysis as in |^ Proposition 7] we 
conclude that the whole procedure requires 0(^(nL 1 + n n )M(D J )Ej n+ i^ arithmetic 
operations in Q. 

Finally, combining the complexity estimates of the three steps above and the 
probability of achievement of the two generic conditions imposed to the coefficients 
we deduce the statement of the proposition. □ 

Using the algorithm of the statement of Proposition 14. 61 for all 7 s T we obtain 
approximations of the factors m 7 which allow us to compute the minimal poly- 
nomial m u and hence a geometric solution of V . Our next result outlines this 
procedure and estimates its complexity and error probability. 

Proposition 4.7. Suppose that we are given a geometric solution of the vari- 
ety ^0,7 for all 7 £ T, as provided by Theorem \4-3\ with a linear form u G 
Q[X\, . . . ,X n ] whose coefficients are randomly chosen in the set {1, . . . , 4pZ? 4 }, 
where p is a fixed positive integer. Then there is an algorithm which computes a 
geometric solution of the curve V with error probability bounded by 1/p performing 
0({n 2 L + n 1+n )M(M r )M(D)(M(D) + M(E))) arithmetic operations in Q. Here 
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L := max 7£ r L 1 , where L 7 is the number of arithmetic operations required to evalu- 
ate the polynomials hi tl of (|3.6(l for a// 7 G T and Mr '■= max 7g r max {71 , . . . ,7 n +i}. 

Proof. For each 7 G T, we apply the algorithm underlying the proof of Proposition 
14.61 in order to obtain an approximation of m u with precision 27„ + i_E. Due to 
(|4.14fl . this polynomial immediately yields an approximation with precision 2E of 
m 7 (T,F) in Q((T 1 /7« +1 ))[y]. 

Multiplying all these approximations, we obtain an approximation with precision 
2E of the polynomial fh u — Il-yer m i °f | ]4.12[l . Since every coefficient a,j(T) of 
fh u G Q(T)[y] is a rational function of Q(T) having a reduced representation with 
numerator and denominator of degree at most E, such a representation of Oj(T) can 
be computed from its approximation with precision 2E using Pade approximation 
with 0(M(E)) arithmetic operations in Q. 

In order to estimate the complexity of the whole procedure, we estimate the 
complexity of its three main steps: 

(i) the computation of the polynomials m 7 with precision 2E for all 7 G T, 
which requires 0( Y, 7 er( nL f + n a )M(D 7 ) (M(M 7 )M(D 7 ) + £?7„ + i)) arith- 
metic operations in Q, 
(ii) the computation of the product Il 7 er m 7 w hh precision 2E, which requires 

0{U(D)M(E)) arithmetic operations in <Q>, 
(Hi) the computation of a reduced representation of all the coefficients of fh u G 
Q(T)[Y], which requires 0(M(E)D) arithmetic operations in Q. 

In conclusion, the algorithm performs 0((nL + n n )M(Mr)M(D)(M(D) + M(E))) 
arithmetic operations in Q, where Mr ■— majCy e r{Mj, 7^+1} and L := max 76 r Ly 
Next we discuss how this procedure can be extended to the computation of a 
geometric solution of V in the sense of Section 1231 Two computations of the above 
procedure involve divisions by the coefficients Ui of the linear form u: the compu- 
tation of the resultant of l|4.15|) for all 7 G T and the Pade approximations of (Hi). 
Both computations are reduced to D applications of the EEA, which is performed 
in a ring Q(A). A similar analysis as in Proposition ^. ll shows that all the denomina- 
tors in Q[A] arising during such application of the EEA are divisors of a polynomial 
of degree 4Z? 4 . Therefore, according to Lemma f2.4l we conclude that a geometric 
solution of V can be computed with O ( (n 2 L + n 1+n ) M ( Mr) M (D) (M (D) + M (E)) ) 
arithmetic operations in Q, with an algorithm with error probability at most 1/ p, 
provided that the coefficients of u are randomly chosen in the set {1, ... , ApD 4 }. □ 

Putting together Theorem 14.31 and Proposition 14 . 71 we obtain the main result of 
this section: 

Theorem 4.8. Let p be a fixed positive integer. Suppose that the coefficients of the 
linear form u of the statement of Theorem \4-3\ and of the linear form u are randomly 
chosen in the set {1, . . . ,4np_D 4 }. Then the algorithm underlying Theorem \4-tl\ and 
Proposition ^. 7| computes a geometric solution of the curve V with error probability 
3/p performing 0((n 2 L + n 1+n )M(M r ) log(Q)M(D)(M(D) + M(E))) arithmetic 
operations in Q. Here L := max 7£ r i 7 , where L 7 is the number of arithmetic 
operations required to evaluate the polynomials hi n of <|3.6[) for all 7 G T, Q := 
2maxi<i< r j{||g||; q G A,;}, and Mr ■= max 7gr ||7||. 
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4.3. Solving a sufficiently generic sparse system. Now we obtain a geometric 
solution of the zero-dimensional variety V\ :— {x € C™ : hi(x) — 0, . . . , h n (x) = 0} 
from a geometric solution of the curve V. 

With notations as in the previous section, we have that V\ = 7r _1 (l), where 
7r : V — * A 1 is the linear projection defined by ir(x, t) := t. Moreover, due to 
Lemma 1431 the equality V\ = 7r _1 (l) n V holds. 

This enables us to easily obtain a geometric solution of V\ from a geometric 
solution of the curve V. Indeed, let m u (T, Y),v\(T, Y), . . . , v n (T, Y) be the poly- 
nomials which form a geometric solution of V associated to a linear form u G <Q[A]. 
Suppose further that the linear form u separates the points of V\ . Making the sub- 
stitution T — 1 , we obtain new polynomials fh u (1, Y) , Vi (1, Y) , . . . , v n ( 1 , Y ) G Q[Y] 
such that fh u (l,u(X)) and ^r-(l,u(X))Xi — Vi(l,u(X)) (1 < i < n) vanish over 
V\. Taking into account that degy(m„) = D = #V± and that u separates the points 
of Vi, it follows that the polynomials m„(l, Y), Vi(l, Y), . . . , v n (l, Y) G Q[Y] form 
a geometric solution of V\ . 

Proposition 4.9. Let p be a fixed positive integer. With assumptions and nota- 
tions as in Theorem \4-8\ the algorithm described above computes a geometric solu- 
tion of the zero- dimensional variety V\ with error probability 4/ p using 0((n 2 L + 
n 1+n )M(A4 r )log(Q)M(I))(M(D) + M{E))) arithmetic operations in Q. 

5. The solution of the original system 

Let notations and assumptions be as in the previous sections. Assume that we are 
given a geometric solution m u (Y), v\(Y), . . . , v n (Y) of the zero-dimensional variety 
V\ defined by the polynomials hi ■= fi + g%, . . . , h n := /„ + g n . Assume further 
that the linear form u of such a geometric solution separates the points of the zero- 
dimensional variety /i = • • • = f n = 0. In this section we describe a procedure for 
computing a geometric solution of the input system f\ = ■•■ = /« = 0. 

For this purpose, we introduce an indeterminate T over Q[A] and consider the 
"deformation" Fi, . . . , F n G Q[X, T] of the polynomials fi, . . . , /„ defined in the 
following way: 

Fi(X,T) := MX) + (1 - T)g t (X) (1 < i < n). 

Set V := {(a?, t) G A" +1 : Fi(x,t) = ■ ■ ■ = F n (x, t) = 0} and denote by vr : V -> A 1 
the projection map defined by n(x,t) := t. As in Subsection 14.21 we introduce the 
variety Vdom C A™ +1 defined as the union of all the irreducible components of V 
whose projection over A 1 is dominant. 

5.1. Solution of the second deformation. In this section we describe an efficient 
procedure for computing a geometric solution of Vdom from the geometric solution 
of 7r -1 (0) provided by Proposition 14. 91 

Since 7r _1 (0) is the variety defined by the "sufficiently generic" sparse system 
hi(X) = Fi(X, 0) = 0, . . . , h n (X) — F n (X, 0) = 0, with similar arguments to those 
leading to the proof of Lemma [4.51 it is not difficult to see that the polynomials 
Fi, . . . , F n , the variety V, the projection it : V — > A 1 , and the fiber 7r _1 (0) satisfy 
all the assumptions of Lemma f4. 41 We conclude that Vdom is a curve and that the 
identity V H 7r _1 (0) = Vdom H 7r _1 (0) holds. Furthermore, Lemma f4 . 41 implies that 
all the hypotheses of (SHI Theorem 2] are satisfied. 
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Therefore, applying the "formal Newton lifting process" underlying the proof 
of [H3 Theorem 2], we compute polynomials m(T, Y), v%(T, Y), . . . , v n (T, Y) <E 
Q[T, Y] which form a geometric solution of Vdom- The formal Newton lifting process 
requires 0((nL' + nP +1 )M(D)M(E'j) arithmetic operations in Q, where L' denotes 
the number of arithmetic operations required to evaluate F\ , . . . , F n and E' is any 
upper bound of the degree of m in the variable T. 

We can apply Lemma l2.3l in order to estimate deg T m in combinatorial terms. In- 
deed, let Qi, . . . , Q n C M" +1 be the Newton polytopes of the polynomials Fx, . . . , F n 
and let A C R n+1 be the standard unitary simplex in the plane {T = 0}. Since 
Qi C Qi x [0, 1] holds for 1 < i < n, where Qi C K. ra is the Newton polytope of hi, 
by H2.3fl of Lemma \l. 31 we deduce the following estimate: 

n 

(5.1) deg T m u <E' :=J2 MV(A, Q u ..., Q^,Q l+1 , Q n ). 

With this definition of E', we have: 

Proposition 5.1. Suppose that we are given a geometric solution of the variety 
V\ , as provided by Proposition \4.iJ[ A geometric solution of Vdom can be deter- 
ministically computed with 0({nL' + n n+1 )M(D)M(E')^J arithmetic operations in 

Q. 

5.2. Solving the input system. Making the substitution T = 1 in the polynomi- 
als m(T,Y),Vi(T,Y) (1 < i < n) which form the geometric solution of Vd om com- 
puted by the algorithm of Proposition l5.1l we obtain polynomials m(l, Y), i>i(l, Y), 
...,v n (l,Y) 6 Q[Y] which represent a complete description of our input sys- 
tem fi(X) = ■■■ = f n (X) = 0, eventually including multiplicities. Such mul- 
tiplicities are represented by multiple factors of m(l,Y), which are also factors 
of i>i(l, Y), . . . , v n (l, Y) (see e.g. |22 §6.5]). In order to remove them, we com- 
pute a(Y) :— gcd(m(l, Y), (dm/dY)(l, Y)) , and the polynomials m(l,Y)/a(Y), 
(dm/dY)(l,Y)/a(Y), Vi(l,Y)/a(Y) (1 < i < n). These polynomials form a geo- 
metric solution of our input system and can be computed with 0(rM{D)E'^ addi- 
tional arithmetic operations in Q. 

Summarizing, we sketch the whole procedure computing a geometric solution of 
the input system fx = ■ ■ ■ = f n = 0. Fix p > 4. We randomly choose the coefficients 
of the polynomials gi,...,g n in the set {1, . . . , ip(nd) 2n+1 + 2pn 2 2 jVl+ " +A ^ S } and 
coefficients of linear forms u, u in the set {1, . . . , 16npD 4 }. By Theorem £21 it 
follows that the polynomials g\,...,g n and the linear forms u,u satisfy all the 
conditions required with probability at least 1 — l/p. Then we apply the algorithms 
underlying Propositions 14.91 and 15.11 in order to obtain a geometric solution of the 
variety Vdom- Finally, we use the procedure above to compute a geometric solution 
of the input system f\ = ■ ■ ■ = f n = 0. This yields the following result: 

Theorem 5.2. The algorithm sketched above computes a geometric solution of the 
input system /i = • • • = f n = with error probability at most 1/p using 

0((n 2 max{i, L'} + n 1+n ) M(D) ( log(Q)M(X r )(M(L>) + M(E)) + M(£'))) 

arithmetic operations in Q. Here L :— max 7e r^7j where L 7 is the number of 
arithmetic operations required to evaluate the polynomials hi n of ()3.6() for allj £ T, 
L' denotes the number of arithmetic operations required to evaluate Fx, . . . , F n , 
Mr ■■= max 7er ||7||, and Q := 2max 1 < 4 <„{||g||; q e AJ. 
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